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ABSTRACT 


I 


A  research  program,  comprising  both  analytical  and 
experimental  tasks,  has  been  pursued  with  the  objective  of 
examining  the  influence  of  biaxial  applied  loads  on  the 
mechanical  behavior  of  bodies  containing  cracks.  / 

Under  the  analytical  part  of  this  program,  £  general 
fracture  mechanics  analysis  was  performed  to  examine  the 
influence  of  biaxial  applied  loads  on  the  mechanical  state 
of  the  body.  The  geometries  examined  were  the  single  crack 
and  two  coplanar  cracks  with  an  arbitrary  orientation  and 
the  cracked  shear  panel.  It  was  found  that  the  biaxial 
loads  influenced  all  aspects  of  the  mechanical  state  of  the 
body,  with  the  exception  of  the  stress  intensity  factor  for 
a  crack  oriented  parallel  to  the  biaxial  load.  The  extent 
and  nature  of  the  biaxial  effect  on  the  crack-tip  stress 
field,  stress  intensity  factor,  angle  of  initial  crack 
extension,  crack-tip  displacements,  elastic  strain  energy, 
fracture  load  and  fatigue  crack  growth  rates  are  all  pre-1'  1 
sented  in  this  report. 

The  experimental  part  of  this  program  involved 
developing  a  biaxial  test  facility’,  and  performing  a  con¬ 
siderable  number  of  photoelastic,  fracture  \toughness  ,  and 
fatigue  crack  growth  rate  experiments.  Confirmation  of 
the  analytical  predictions  was  obtained  for  the  biaxial 
effects  on  the  crack-tip  stress  field,  the  angle  of  initial 
crack  extension,  the  fracture  load,  and  the  fatigue  crack 
growth  rates.  The  biaxial  loads  were  seen  to  influence 
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SECTION  1 


INTRODUCTION 

Ever  since  the  seminal  research  of  Griffith  into  the 
nature  of  the  fracture  of  solids  in  early  1920,  the  opinion 
of  many  associated  with  the  materials  science  and  fracture 
mechanics  communities  has  been  that  loads  applied  parallel  to 
an  existing  crack  (in  addition  to  the  tensile  load  perpendi¬ 
cular  to  it)  have  no  effect  on  fracture  behavior.  This  idea 
was  first  advanced  by  Griffith  himself  in  postulating  his 
well  known  criterion  for  the  onset  of  crack  instability  [12, 

13],  The  mathematical  or  analytical  justification  for  such 
an  assertion  can  be  shown  to  be  lacking  however,  and  has  been 
a  source  of  controversy  over  the  years  [48], 

Apparent  corroboration  of  this  idea  of  the  unimportance 
of  loads  applied  parallel  to  a  crack  (in  a  biaxial  load  situ¬ 
ation)  appeared  to  come  some  thirty-five  years  later,  when 
Irwin  introduced  the  crack-tip  elastic  stress  intensity  fac¬ 
tor  approach  to  fracture  mechanics  [49,50].  Irwin's  deriva¬ 
tion  of  the  singular  one-term  series  approximations  for  the 
elastic  stress  in  the  neighborhood  of  the  crack-tip,  utilized 
Westergaard ' s  form  for  the  solution  of  the  problem  of  a  biax- 
ially  loaded  plane  infinite  sheet  with  a  central  horizontal 
crack  [51],  Westergaard ' s  solution,  however,  incorrectly 
omitted  a  term  containing  the  horizontal  load  parameter,  (i.e., 
the  second  term  of  the  right  side  of  Eq.  (2.11)  of  this  report). 
Consequently,  in  Irwin's  subsequent  series  representations  of 


the  stress  field  near  the  crack-tip,  no  manifestation  of  the 
presence  of  the  horizontal  load  could  possibly  appear.  How¬ 
ever,  Irwin  seemed  to  understand  the  significance  or  importance 
of  loads  parallel  to  the  crack,  because  he  added  a  constant 
but  unspecified  term,  ctqx>  to  his  expression  for  the  x-direction 
stress  component  near  the  crack-tip,  [SO].  It  is  now  under¬ 
stood  that  this  term  corresponds  to  the  second,  non-singular, 
term  of  Eq.  (2.24)^  of  this  report. 

Thus  two  analytical  errors,  with  the  latter  one  accidently 
having  the  effect  of  masking  the  presence  of  the  earlier  one, 
led  to  a  general  misunderstanding  of  the  significance  of  the 
biaxial  load  effect.  This  situation  still  persists.  Now 
that  fragments  of  contradictory  experimental  data  have  begun 
to  appear,  during  the  last  five  to  ten  years,  the  need  to 
clarify  this  question  is  apparent. 

This  research  effort  addresses  itself  to  this  problem, 
that  is,  to  the  question  of  the  presence,  or  lack  of  presence, 
of  the  effect  of  load  biaxiality  on  fracture  or  on  fracture 
related  behavior.  It  attempts  to  examine  the  problem 
systematically  and  in  detail,  utilizing  specially  designed 
experiments  where  possible  to  check  on  the  validity  of  the 
analytical  findings.  The  experiments  were  performed  on 
materials  that  are  used  extensively  in  aircraft  structures. 


) 


-2- 


PART  I 


ANALYSIS 

When  a  cracked  body  is  subjected  to  forces  that  act 
parallel  to  the  plane  of  the  crack  fin  addition  to  the  ten¬ 
sile  forces  applied  perpendicular  to  it),  any  influence  that 
the  parallel  loads  may  have  on  the  fracture  behavior  of  the 
body  should,  if  present,  manifest  itself  in  several  different 
ways.  Moreover,  if  the  biaxial  load  effect  exists  and  is 
significant,  it  must  obviously  affect  the  mechanical  state  of 
that  region  of  the  body  in  close  proximity  to  the  crack  tip. 
Investigation  of  the  problem  should  therefore  direct  initial 
emphasis  on  illuminating  the  possible  relation  of  load  bi- 
axiality  on  those  aspects  of  fracture  that  can  be  observed 
experimentally.  These  would  include  crack-tip  shear  iso¬ 
chromatics,  which  may  be  seen  photoelasti cal ly ,  the  angle  of 
initial  crack  extension,  which  may  be  observed  directly,  as 
can  the  rate  at  which  fatigue  crick  growth  takes  place  under 
cyclic  normal  loads  at  varying  states  of  load  biaxiality. 

The  critical  value  of  the  normal  tensile  load,  that  is,  the 
fracture  load,  can  also  be  measured  at  different  levels  of 
load  biaxiality,  which  suggests  that  the  analysis  should 
concern  itself  with  this  facet  of  the  problem  as  well. 

The  infinite  sheet  geometry  was  chosen  for  the  theoret¬ 
ical  part  of  the  investigation  because  exact  solutions  for 
several  cracked-body  geometries  are  known,  thereby  permitting 
derivation  of  specific  results  that  give  explicit  demonstration 
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of  the  effects  of  the  load  applied  parallel  to  the  plane  of 
the  crack.  Although  experiments  must  be  performed  on  test 
specimens  that  have  finite  dimensions,  analysis  based  upon 
infinite  body  geometry  can  be  useful,  nevertheless,  because 
it  can  serve  as  a  prototype  that  provides,  at  the  very  least 
valuable  information  of  an  insightful  and  qualitative  nature 


SECTION  2 


STRESS  AND  DISPLACEMENT  IN  THE  VICINITY  OF  THE  CRACK  TIP 

In  the  linear  theory  of  isotropic  elastic  solids  in 
which  body  forces  arc  neglected,  the  plane  problem  of  statics 
is  reduced  to  a  boundary  value  problem  for  the  biharmonic 
stress  function.  The  same  problem  can  also  be  formulated 
equivalently  in  terms  of  a  pair  of  analytic  functions,  $  and 
ft ,  of  the  complex  variable  z»x+iyBre “  ,  which  uniquely  specify 
the  plane  stress  and  displacement  components  by  means  of  the 
following  relations  [1]: 

txx  +  tyy  ■  2C<trz)  +  *  ( 2 )  1 

tyy  -  txx  *  2itxy  -  2((z-z)||  ♦  no  *  *('z)l  (2.1) 

2w[ux  ♦  iuy]  «  «P(z)  -  w(z)  -  (z-z)TCzl  , 
where  the  ove"b..r  designates  complex  conjugation  and 

4>(  z)  «y*<t(r)dz,  «(z)  *  f  ft(z)dz  .  (2.2) 

When  for  a  given  boundary  value  problem,  that  is,  for  a  given 
applied  load  anc  cracked  body  geometry,  the  functions  $  and 
ft  are  determined,  then  the  stress  and  displacements,  and  all 
other  mechanical  quantities  expressed  in  terms  of  them,  are 
determinable  from  Eqs.  (2.1)  and  (2.2). 

2.1  SINGLE  CRACK 

The  plane  components  of  the  stress  tensor  referred  to 
the  pair  of  coordinate  systems  shown  in  Fig.  1  are  related 
by  expressions 
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t  -  t  ♦  2it  •  e2i3(t'  -  t'  *  2it '  ) 

yy  xx  xy  1  yy  xx  xy' 


•  “<*«  ‘  ^y  *  2itiy> 


(2.3) 


t  +  t  *  t  +  t 
yy  xx  yy  xx 


(2.4) 


The  boundary  condition  along  the  exterior  boundary  of  the 
plane  infinite  body  having  an  inclined  interior  crack  geometry 
(cf.  Fig.  1)  is  specified  to  be 


*  cr  , 


t;x(«)  -  ko  ,  t-y(»)  *  0  ,  (2.5) 


which  can  be  expressed  relative  to  the  (x-y)  coordinate 
system  by  means  of  Eqs .  (2.3)  and  (2.4)  as 

tyy(-)  »  f(l+k)  -  §(l-k)cos2a 


txx(«)  *  f(l+k)  +  f(l-k)cos2o 


(2.6) 


txy(°°)  »  -^(l-k)  sin2a  , 

Along  the  interior  (crack)  boundary 
tyy(x,0)  »  txy(x,0)  -  0  , 


x  <  a  . 


(2.7) 


The  solution  to  the  boundary  value  problem  posed  by  Eqs. 
(2 . 5) -  (2 . 7)  has  the  form  [2] 


4(z) 


Cl(z) 


.  1X2MBU  .1, 
2  U2-a’)h  2 


where 


|(2A+B)  -  |  (l-e2lot)  +  k(l  +  e2l0t) 


■  -j(l-k)e- 


(2.8) 


(2.9) 


! 


Sir 


From  Eqs.  (2.2),  (2.8)  and  (2.9),  it  follows  that 

♦  (*) 
w(z) 


»  |  (l-e2lct)+k(l>e2l!jt)j  (22-a2)‘s*(l-k)e2ia: 


(2.10) 


When  the  crack  has  a  horizontal  orientation,  a*fr/2,  the 
above  expressions  reduce  to 
Hz)  ' 


ft(8) 

and 

♦(z) 

«(2) 


02 


-  a 


(l-k) 


(2.11) 


If^-a2)1*  ♦  j(l-k) z  . 


(2.12) 


Expression  of  the  stress  and  displacement  components  in 
the  immediate  vicinity  of  the  crack  tip  is  facilitated  by  the 
following  coordinate  transformation  (cf.  Fig.  2) 

z  -  a  »  c  -  re19  ,  (2.13) 

relative  to  which  <t(z)-*-<t(c)  may  be  shown  to  have  the  power 
series  representation 


♦  CO 


in  which 


a(ir)h 

c 

(2^ 

f 

c 

(2ttO^ 

4 

(l-e2i“) 

(f)  *  ^  n-k)e2iaQ°*  |(f J  - 


(2.14) 


2  iou 


However 
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(2. IS) 


.  1,  0  (^3  (1+k)  -  (l-k)cos2o-i  (l-k)5in2a 


where 


KriK2  "  iia[2(27r),S(2‘a)!s^(2)]  C2,16) 

is  the  complex  elastic  stress  intensity  factor  [3,4] .  The 
equivalence  of  expressions  (2. IS)  and  (2.16)  can  be  shown  by 
direct  calculation  of  the  limit  Eq.  (2.16)  with  4  as  given  by 
Eq .  (2.8).  Thus 

h 


K,-iK0  , 

HO  «  — - ♦  j(l *k) e2lot  4.  0(£ 

2(2trO^  4 

and  for  0<  |  q/sl  |  -  (r/a)  «1 ,  we  have  the  following  series 
approximations  for  4  and  taken  to  order  (r/a)*5: 

24(0' 


(2.17) 


2SUO 


..  K.-iK7  -i  4 
■  1  ^  a  i  f(l-k)e 


2ia 


t2;;:io  (t)1  *  K*)3/: 


(2irr) 

Similarly,  a  series  expansion  for  4(0  has  the  form 
•«>  •  ^1©'  • 

•  •  ¥<*■«(>*£)•“*•  °©‘” 
which  can  be  approximated  by 

^  i(lL  -  ie 


(2.18) 


(2.19) 


4(0 

w(0  J 


(K 


j-iK, ){^y  e'v2,J  ^Cl-k)'a*re 


2  iot 


(2.20) 


when  0<|c/a|  ■  (r/a)<<l. 


The  stress  and  displacement  components  near  the  crack  tip 
are  thus  obtainable  from  Eqs.  (2.1),  (2.18)  and  (2.20),  where 
in  Eq.  (2.1)  z  and  z  are  replaced  by  t  and  Q,  with  the  results 


xx 


cos  |  .inf (j-cfco.^) 


♦  o(l-k)cos2a 


K1  9 

'  7^7* cos  1  \ 


( l  +  sin-|sin-y-^  + 


2  !jsin|cos|cos^  (2. 21) 


K 


tjcy  (2ttt) 


1  .9  6  39 

TJ  Siny  cos7  COS-J- 


(2irr) 


(2irr) 


2  ^  cos|-  ^1  -sinysin^p 


) 


and 


u 


cosf[l(K'1)','sin2!]^(H)'i  sinf^(^I)*cos2f 


x  '  TflTf)  “3 

,  U-k)c  r  (cos  (9+ 2a)  +kco:>  (6-2ot)  -  2sin8sir.2cx]  +  (K.  +  l)a  cos 2a 
oU 


(2.22) 


uy  *  t(hV  sinf[|«*l)-cos^]*  cosf[f(l-K)**in2f] 

*d*k)a  I  r  [sin(2a-9)+<sin(2a+8)-2sin0cos2a]  +  (ic+l)a  sin2a|  . 
8U  | 

In  the  above  expressions  the  elastic  stress  intensity  factor 
Xx  and  K2  have,  from  Eq.  (2.15),  the  explicit  forms 

Kx  .  2ll|ll^i+k)-(l-k)cos2aJ 


- ^  -  (l-k)sin2a 


(2.23) 
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The  variation  of  Kx  and  K2  with  the  angle  of  crack  inclination, 
a,  and  the  biaxial  load  factor,  k,  are  shown  in  Figures!  3  and 


For  the  crack  with  a  horizontal  orientation,  a=rr/2, 
Eqs.  (2.21)- (2.23)  become 


(2irr) 


cos||l-sin|-sin~j  -  (l-k)o 


(2rrr 


^  C°sf(l*3in§Sinf) 


(2.24) 


txy  5  - sin§  cos|  cosf- 

xy  (2irr)  J  L  2  2 


uy  =  co.H]*IH34^Sla 


r  cos 0  +  a] 


r  sine 


where 


(2.25) 


h  *  *  K2  -  0  .  (2.26) 

From  examination  of  the  power  series  representations 
given  by  Eqs.  (2.14)  and  (2.19)  for  $(0  and  d(r),  respectively 
it  is  observed  that  the  biaxial  load  factor,  k,  appears  as  a 
coefficient  in  only  the  second  term  of  the  expansions.  Con¬ 
sequently,  in  order  to  take  proper  account  of  the  presence  of 
the  horizontal  load  in  the  approximate  expressions  for  the 
elastic  crack-tip  stress  and  displacements,  a  two-term 


approximation  of  the  series  expansions  for  and  V  must  be 
employed,  as  shown  by  Eqs.  (2,17)  or  (2.18)  and  (2.19)  or 
(2.20).  This  gives  rise  to  the  a (l-k)cos2a  contribution  to 
Eqs.  (2. 21)1  and  (2124) lf  and  to  the  second  lines  of  the  Eqs. 
(2.22)  and  (2.25),  contributions  which  do  not  appear  in  the 
usual  expression  of  these  equations  in  the  fracture  mechanics 
literature.  When  the  crack  is  inclined  the  horizontal  load 
will  also  effect  Kj  and  K2 ,  i.e. ,  will  also  appear  in  the  co¬ 
efficient  of  the  first  term  of  the  series  expansions,  while 
for  phe  horizontal  crack  orientation  is  independent  of  ka 
and  K2at0. 

It  is  also  noted  that  the  customary  expressions  for  the 
displacements  (which  consist  of  only  the  first  lines  of  Eqs. 
(2.22)  and  (2 . 25)) predict  zero  displacements  for  the  crack 
ends  which,  from  physical  considerations  alone,  must  be  in¬ 
correct.  According  to  Eq.  (2.22) 


1-k)  (<+l) aa  (cos  2a 
8 u  /sin  2a 


(2.27) 


whereas  for  the  crack  with  a  horizontal  orientation 


’1-k) (K+l)oa 


r»0 

a-Tr/2 


>  (O 


(2.28) 


a®  tt /  2 


showing  that  the  crack  tips  move  horizontally  out  (away  from 
the  center)  for  horizontal  tensile  loads,  k>l ,  and  move  hori¬ 
zontally  in  for  tensile  loads,  0<k<l,  and  compressive  loads, 
k<0 . 


The  equations  obtained  above  show  explicitly  how  in 
situations  of  biaxial  load  application  along  the  exterior 
boundary  of  the  body,  the  usual  approximation  for  the  stress 
and  displacements  in  the  vicinity  of  the  crack  tip  by  means 
of  only  the  first  term  of  the  series  expansions  of  the  func¬ 
tions  representing  them  is,  in  general,  inadequate ,  because 
it  excludes  in  part,  if  not  entirely*,  the  presence  of  the 
horizontal  load  kc  in  its  representation  of  the  mechanical 
state  of  the  crack  tip  region. 


Referred  to  the  polar  coordinates  at  the  crack  tip  (cf. 
Fig.  2)  the  plane  elastic  strain  energy’ density  (per  unit 
volume)  at  any  point  can  be  determined  by, 


Mr  ,0) 


2  l  xx  3x  V  3y  'xy\,Ty 


+ 


1  f9Ux  cosQ-3  Ux  sin6)  -  t  sine+^i-  cose 

2  [^xATr  COs8  Te  T“7  tyy\~9r  sln0  Te  T~ 

*  r  (9ux  fli,,?ux  cose  x  9 uy  3uy  sineYI 
txy\TTr  5ine*~W  T-" Ir  0  0  50  Te  Tr“/J . 


C2.29) 


By  means  of  Eqs.  (2.21),  (2.22)  and  (2.29)  the  strain  energy 
density  in  the  immediate  vicinity  of  the  crack  tip  is  given 


by  [2,5] 


*  For  the  horizontally  oriented  crack  the  presence  of  the 
horizontal  load  kc  appears  only  through  the  second  term 
of  the  series  expansions. 
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stu?  [*2sin  7  +  (3'<}  sin2  7  +  (K‘1). 


k  ^  r  * 

•g—-.-  I  8sin7|-8sin6'|‘10si,n5|-+16sin4-|+‘2sin2|-+  f<-9>sin2|'+2 

*  .  .  J 


k  k  r  (  / 

+  yg-',"-|  cosf  r64sin9|+128sin7|+(8ic- 25)  sinS|-6  (<+5)  sin3|-+2  (9- <) sin 


+  1 28s in^y-  256s in7|-+l 5  2s in^-|-  24sin^ 


in3f] 


Ki(l-k)a  r  0  ( 

j^t?[C05A 


(tc-l)  *  6sin2|-+8sin^|' I  cos  ;  2a 


K  2  ( 1  -  k )  cr 
4y  (Zttt)'1 


+  2sin||4(<-'2)sin4|+(ll-S(c)sin2|+(K-3)jsin  2a 

-tj  sin-y  -  (<  +  3) -8sin44+l9sin'^'|  j  cos  2a 
)  e  ■  / 


2'cosy  <2  (<-3)  sin4y+  ( 1  -  3  k  )  s  i  n ( ic  + 1 )  T  sin  2a 


a-k)3  f , 
hsir*  [ f 


(  k  + 1 )  c  o  s  -  2  a  -  s  i  n  2  0  s  i  n  4  a 


:4a  j  . 


(2.30) 


Were  one  to  calculate  'p  from  a  quadratic  representation 
entirely  in  terms  of  the  plane  stress  components,  i.e. , 


1  ct2 


1+v  v 


<p  «  *  (t^  +  t  )  +  -V2.  t*  -  ~  t  t 
v  2E  xx  yyJ  E  xy  E  xx  yy 

and  use  Eqs,  (2.21)  with  the  o(l-k)cos2a  term  omitted,  then 

the  resulting  expression  for  \\>  would  have  the  form 

i|/  *  (6)  *^2^2^ 2  (0)  +ci 2K1K2F12  ’ 


-13- 


I  . 


which  is  the  expression  used  in  Sih's  strain  energy  density 
theory  [6).  A  comparison  of  this  approximation  with  Hq.  (2.30) 
shows  that  the  last  three  terms  of  (2.30),  which  stem  from 
the  presence  of  the  non-singular  contribution  a(l-k)cos  2a 
to  Eqs.  (2.21)  and  the  second  line  contributions  to  the 
displacement  equations  (2.22),  are  necessary  additional  terms 
which  cannot  be  ignored  if  proper  account  is  to  be  taken  of 
the  effect  of  the  horizontal  load  ka  on  the  local  (crack  tip) 
elastic  strain  energy  density. 

For  the  horizontal  crack  Eq .  (2.30)  reduces  to 

(*)a/rr2  S  fpf  |_-2sin4f*C3-iOsiii2f*('<-l) 


'  '  \  /  ,  ' 


(1  -  k)'|^y  cosj |(<-l)  +  8sin4y-6sin2j 


3m 


(2.31) 

The  local  elastic  strain  energy  rate  for  an  arbitrarily 
small  circular  region  of  unit  thickness  centered  at  the  crack 
tip  with  radius  Tq,  where  0<rQ/a<<l.  is  obtained  as  the 
derivative  with  respect  to  the  crack  size  of  the  integral 


nr 


^(r,9)rdi'd8  , 


which,  upon: utilisation  of  Eqs.  (3,30)  arid  (2.23),  yields 

•  i  -  - 

~|(2k-1)  J  (1+k)  •  (l-k)cos2a  j  +2(02)  (1  -k)  2sin22c<j 


3a 


d^Ttr, 


32 


[ 


3 (5k- 7)  (1+k) - (1 -k)cos2 


■  14  - 


aj  cos  2 


(2.32) 


a+(6<+l)  (l-k)sin2a }  ,, 


!]. 


I 


. . 


m.-  nmmm  rf.. .  i  iw  ■«<«>  '■  *•<  ■ 


•Vi 


Ijr?-': 


( 


The  variation  c£  (3473a)  with  crack  angle  and  load  biaxiality 
is  shown  in  Fig.  S  for  a  value  of  u«0.30  and  (r0/a)“0.Q7. 

For  the  horizontal  crack  Eq.  (2.32)  simplifies  to 

h 


a®  it/  2 


C2TTr( 

1 6  u 


(2k- 1) 


o  i6js^7iLi.:.k2 

1 2a  /  1  Sir 


(2.33) 


2.2  TWO  COLLINEAR  CRACKS 

For  the  double  collinear  crack  geometry  shown  by  Fig.  6, 
the  exterior  boundary  conditions  are  specified  by  Eqs.  (2.6), 
whereas  along  the  interior  crack  boundaries 

t y y ( x , 0 ) ™ t Xy ( x , 0 ) *  0,  | a | < | x | < | b |  .  (2.34) 

The  solution  is  given  by  the  sectionally  holomorphic 
functions  [7] 

*  ( z ) 


G(z) 

♦  (z) 
w(z) 


_ Mfr. M  ..  ,  .  c 

((z2-a2)(z2-b2)]h 


(2.35) 


■/ 


_ z2dz  _ 

[(z2-a1)(zZ-bZ))li 


-  AB 


/ 


_ dz _  .  rT 

[(zz-az)(z2-b2))h 

(2.36) 


where 


*[ 


(l-e2la)+k(l+e2ia) 


> 


B  *  b2[E(p)/F(p)]  -  b2 A (p) ,  p2»l-a2/b2 

C  -  'jo(l -k) e2la  . 


(2.37) 
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Note  that  in  the  limit  as  a*0  the  complete  elliptic  integrals 
E (p ) -►  1  *  F(p)-*"»  and  the  above  expressions  reduce  to  the  solution 
(2.8)  and  (2.10)  for  a  single  crack  of  length  2b. 

In  establishing  the  structure  of  the  functions  $  and  ft 
near  the  ends  of  the  cuts,  because  of  the  geometrical  symmetry, 
only  one  cut  need  be  examined,  say  the  one  along  the  positive 
x-axis.  Consider  first  a  circular  region  enclosing  the 
end  of  the  crack  at  z*b,  with  radius  sufficiently  small  so 
as  to  exclude  the  point  z«a,  (cf.  Fig.  7).  Equation  (2.35) 
can  be  rearranged  to  the  form 


Hz) 
fl(  z) 


1 

(z-b)Fi 


r  az2 

^(z+b^tz^a2)4 


AB 

(z+b^fz^aV* 


C 


fb(z) 

cV-b? 


*  C, 


where  f^(z),  representing  the  terms  appearing  within  the 
bracket,  is  clearly  analytic  for  all  points  z  in  Rb .  A 
similar  argument  can  be  made  for  all  points  of  a  circular 
region  R  enclosing  the  end  of  the  cut  at  z*a,  while  excluding 

the  point  at  z«b,  merely  by  interchanging  b  with  a  in  the 
above  expression.  Designating  the  ends  of  the  cut  by  z»z^ * 
j*l,2,  so  that  z^*a  and  Z2*b,  the  above  arguments  can  be 
jointly  stated  as 


Hz) 

Q(z) 


H  (z) 

*  c’ 


(Z'Zj  ) 


j*l  ,2 
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where  the  functions  fj (2)  are  analytic  for  all  points  2  of 
the  regions  Rj ,  thus  allowing  the  series  representations 

f.  (z)  -  sj  +  sj(Z-Zj)  r  0^  (2- )  2  *  . 

in  which  the  coefficients  8“[  include  the  constants  A  and  B. 

n 

Thus  in  the  vicinity  of  the  end  points  z^  of  the  cracks  the 
structure  of  and  fi  have  the  form 


G(Z) 


8 


( 2  ■  2  j  )  ^ 


*  3J1(2-2j)li  ♦  8J2(2-2j) 


3/2  + 


*  C 


( 2  ■  2  j  ) 


4  *  c  ♦  0(  |  z-z.  |*)  , 


ZjCRj,  j*l,2.  (2.38) 


Equations  for  the  elastic  stress  and  displacements  in 

“he  crack-tip  region  can  now  be  expressed  in  the  standard 

form.  This  may  be  accomplished  most  economically  by  intro- 

dt  tion  of  a  pair  of  polar  coordinate  systems  originating  at 

i  0 

'  ends  of  the  crack  by  writing  z-Zj=re  ,  with  the  stipulation 
that,  at  the  end  Z2*b  the  angle  0  is  to  be  measured  counter 
clockwise  from  the  real  axis,  whereas  at  the  end  z^*a  it  is  to 
be  taken  clockwise.  Accordingly,  from  Eqs.  (2.1)  and  (2.38) 


t  +  t 
yy  xx 


4Re  C  +  O(r’) 


■yy 


t  xx  +  2ltxy 


-i  se 

-g  e  2  sine 


(2.39) 


+  4Re  C  +  0  ( r*5 ) . 

Consistent  with  standard  notation,  the  coefficient  3q  can  be 
expressed  in  the  form  of  a  complex  elastic  stress  intensity 
factor  as 
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3i  -  - — r  (K,  .  -  IK,.  )  . 

0  2  (270*  V  lj  2j' 


(2.40) 


Equati  as  (2.39)  and  (2.40)  then  yield  the  characteristic  forms 


:  •'  -X- 

xx  Tw* 


K, 


cos4^1'Sin|-sin^r^-  - siny ^2+cosyCOS^^ 


V  ‘  ^  eos^l.slnfsin^J.  sinf  cosf  cos^  .  (2.41) 


4_e  .  3e 


) 


♦  a(l-k)cos2a  . 

K2j  ,4-1  —  —,36 


K,  . 


(2ttt  ) 

K,  . 


xy 


( 2-rrr ) 


^  sinf  C04  C0ST  +  COSt  (1*sinf  sinx) 


( 2  rrr ) 


for  zeRj  ,  j  *1 , 2 . 


At  the  inner  and  outer  ends  of  the  cut,  respectively, 

Klj  *  K11  5  Ki^5  •  *ij  ■  t12  2  *i(l» 

K9j  a  =  K2  (a)  ,  K  2  j  *  ^2">  =  k^fk)  .  (2.42) 

Considering  the  manner  in  which  angle  9  is  defined,  Eqs.  (2.38) 
and  (2.40)  afford  explicit  determination  of  the  stress  intensity 
factors  according  to  the  relations 


(v11^ 


lim 

z+z 


z  ^2  (27r)*(z-Zj  )*t(z)  j  , 


(2.42a) 


from  which 


K.(a)  -  AW'a>l  [  (1  +  k)  -  (1  -k)  cos2a] 

1  2a(b  -a  j”* 


ot^Tiafb^A  (p)  -a2] 

K2(a)  *  2a(b2-a2)J‘  (l-k)sin2a 


(2.43) 
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Kx  Cb) 


K2  Cb) 


-  g^rb  bfl-AC] 
2(b*-aT 

,  £  b  (1-  A  C| 
2 (b  -a  j 


[  (1+k) - (1 -k)cos2a] 


(1 -k) s in2ot 


(2.44) 


for  0<a<b.  The  variation  of  K^(a)/K^(b)  and  K 2  (a ) / K 2  Cb )  with 
crack  spacing  is  shown  in  Fig.  8.  When  the  distance  between 
the  cracks  is  the  equivalent  of  the  crack  length,  interaction 
of  the  inner  crack-tip  stress  fields  is  just  about  negligible 
K1(a)  and  K2(a)  being  only  two  percent  greater  than  K1(b)  and 
K2(b).  For  crack  spacing  at  about  three  tenths  the  crack 
length  or  less,  the  stress  intensity  at  the  inner  crack  ends 
increase  rapidly  relative  to  the  outer  ends. 

For  horizontally  oriented  cracks. 


Kx(a) 


Kx(b) 


a  ( b  -  a  ^  j* 

a  (itb)  ^b  [  1  -  A  (t 
(b2 -a)^5 


(2.45) 


K2(a)  *  K2(b)  -  0 

independent  of  the  load  biaxiality. 

From  Eqs.  (2.1),  (2.2)  and  (2.38),  $(z)  and  w(z)  can  be 
represented  by  the  series 


<Hz) 


oj(z) 


23j(z-zj)1*  *  Cz  +  0(|z-Zj  |3/2) 


(2.46) 


Substitution  of  Eqs.  (2.46),  (2.38)  and  (2.40)  into  (2.1)^ 
yields  the  following  expressions  for  the  displacements  com¬ 
ponents  in  the  immediate  vicinity  of  the  ends  of  the  cracks: 

ux  *  -^(h)  C05f[l(,c'1)*sin2l]*-n1(?r)  sinf[^tx*l)»cos2|] 


-j(<+1)  ZjCOs2a-rsin0sin2a+jr  }<cos  (e+2a) +cos  (0- 2a)J 


(2.47) 


uy  5  -^(h)  sin|[|(xn)-cos2|]  ^(h)  cos|[^a-Kl*5in2f] 

*  ^(*+1) 2  j  sin2a-  rsin9co52a»yr  |>csin(2oc»-9)+sin(2ot-8)jj  , 

where  j»l,2  and  z^a,  z?«b. 

The  same  comments  that  were  made  for  the  single  crack 
concerning  the  presence  of  the  horizontal  load  ko  in  the 
approximations  for  the  crack  tip  region  stress  and  displace¬ 
ments,  apply  verbatim  for  the  double  crack  geometry,  as  Eqs. 
(2.37),  (2.38),  (2 .41)- (2.47)  indicate. 


2.3  SHEAR  PANEL 

The  centrally  cracked  sheet  subject  to  pure  shear  along 


the  outer  boundary  (cf.  Fig.  9) 

tyy(°°)*0»  txx^°°)*0,  txyf<”-)"T 


(2.48) 


relative  to  the  (x'~  y“)  coordinate  system,  can  be  transformed 
to  a  biaxial  boundary  condition 

tyy  (°°)  *-xsin2  8,  txx(»)*Tsin2B,  t  (•) -tcos2  0  (2.49) 

relative  to  the  (x-y)  coordinate  system.  The  boundary  con¬ 
dition  along  the  interior  crack  boundary  is  specified  by  Eq.  (2.7)  ) 


I 


l,„-2is/  z  f  A 
71Te  l  TTZ'TJT  l)  * 


The  solution  to  the  above  problem  is  given  by  [8] 

♦  (2)  ' 

r  i  .^a/  ■» 

*  1 

a(z) 

so  that 

♦  (z) 

(  -  _14 _ -2ifl 

Ll 


(2.50) 


“>(*) 


■] 


(z  -a“) 

Relative  to  the  polar  coordinate  system  at  the  crack  tip 

*ies  expansio 

M 


(2.51) 


z-a*s*rel6  (cf. 

Fig.  2), 

4>  has  the  power 

$ r r}  *  -iire'2^^ 

1  1/iY 

”  •  <«'(* )°  * 

J  *  -y1Te 

K.  -iK, 

=  — - r  + 

|ixe'2iS 

*  o(|f I*1)  , 

2  (2-rrc)  ^ 

\  a  / 

*la 


(2.52) 


where  by  means  of  Eqs.  (2.16)  and  (2.50) 

K1-iK2  *  -ix(7ra)V2i8  .  (2.53) 

Thus  for  0< | s/a  I  * (r/a) <<1 ,  up  to  terms  of  order  (r/a)  ,  the 
stress  functions  become,  (analogous  to  Eqs.  (2.18)), 

2$  (S) 


2H(S)J 


Kl~lK2  ^-i (0/2)  . 
(2ttt)^ 


±  ire 


2  i  0 


Similarly,  for  <J>(s)  and  u(s)  up  to  terms  of  order  (r/a) 
(analogous  to  Eqs,  (2.20)), 


(2.54) 

3/2 


-2iS 


(2.55) 


The  elastic  stress  and  displacement  in  the  vicinity  of 
the  crack  tip  follow  from  Eqs.  (2.1),  (2.54)  and  (2.55)  as 
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\ 


t  * 

yy 


txy  a 


( 2  ttt 


(2irr) 


(2irr) 


7^  cos?( 


1-sinf.in^)-— ^  sinf( 


2+cos|cos:y-  1  +2Tsin8 


cos^ll  +  sin-ysin-y^*- - — ^  siny  cos|*  cos4p 


siny  cos  !>  cos-^p  +  - r  cosy 

7  2  “  (2  ITT)*5  7 


-  si 


iny  sin4p^ 


(2.56) 


s  t(h)  cosf[T(<‘l)*sin2!]*lf(T?)  sinf[T<,<*1’*cos2|] 

*  Ju^ CK+^) (rcos8+a) sin2  8- (k+3) rsin9cos2  3] 

S  t(&)  sinl[TfK*1)-<:os2!]’lf(Ti-)  coS|[|(l-,c)*si„2|] 

+  ~j[ (<+l) (rcos0+a) cos2 3+ (<-5) rsinQs in2 3]  , 


where  from  Eq.  (2. S3) 


■  -t  (-rra)  1*sin2  3  ,  (Tra)>scos2  3 


(2. 57) 


(2.58) 


It  will  be  instructive  to  again  focus  attention  on  the 
second  term  of  the  series  expansion  approximations  given  by 
Eqs.  (2.54)  and  (2.55),  which  lead  to  the  2tsin3  contribution 
to  Eqs.  (2.56)  for  the  stress,  and  to  the  second  lines  of  Eqs. 
(2.57)  for  the  displacement.  As  for  the  biaxially  loaded 
cracked  bodies,  failure  to  retain  these  terms  in  the  series 
approximations  for  the  crack  tip  stress  and  displacement  can 
lead  to  significant  error  in  describing  the  mechanical  state 
about  the  crack  tip.  For  example,  for  the  displacements  at 
the  tip  of  the  inclined  crack 
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u  J 

f  2  TC<n)a 

j  sin2 3 

u 

j  cos2B 

• 

o 

which  according  to  the  usual  one-term 
should  be  zero. 


» 


series 


(2.59) 

approximation 
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SECTION  3 


ISOSTATS  FOR  MAXIMUM  SHEAR  STRESS 


The  maximum  shear  stress  at  any  point  of  a  body  in  plane 
stress  or  plane  strain  is  determined  by 


2  _  1  /_  .  n2  _  .2 


m 


4^yy  txx 


y  +  t 


xy 


(3.1) 


3.1  SINGLE  CRACK 

For  the  maximum  shear  stress  near  the  tip  of  the  inclined 
crack  Eq.  (2.21)  may  be  substituted  into  Eq.  (3.1)  with  the' 
result 


.2  s 
m 


i^^K12sin20  +  K22  (4-3sin29)+4K1K2sin9cos|-J 
ya(l-k)cos2aj" 


L(2irr)H 

:.]  . 


3  0 

sinQs 


K2 


(Sl 


3  6 

in0cos-y-+2s 


inl) 


ya(l-k)cos2a  |  .  ■  (5.2) 

With  the  crack  horizontally  oriented,  the  maximum  shear, 
stress  becomes 


(X-L, 


jiyK12sinz8+- 


K, 


-sin0sin-^o(l-k)+j  (1-k)2.  (3.3) 


r/2  °"A  A  2  (2irr) 

If  the  customary  singular  approximation  for  stress  were 
employed  in  Eq,  (3.1),  i.e.;,  Eq.  (2.21)  without  the  o(l-k)cos2a 
contribution,  only  the  first:  line  of  Eq.  (3.2)  would  appear 
for  the  inclined  crack  and  only  the  first  term  of  (3.3)  would 
appear  for  the  horizontal  crack.  Designating  the  latter  by 


(T"l 


tt/ 2 


8W  Kl2  sinh  ’ 


(3.4) 
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2  2 

where  a  ira  *  it  is  seen  that  the  effect  of  the  horizontal 
load  ko  on  the  maximum  shear  stress  disappears  entirely.  The 
ratio,  of  Eq.  (3.3)  to  Eq.  (3.4) 


&) 


a*7r/  2 


Msin^-M  (l-k^-^yj  sin9sin-^-+2  (1-k) 

(f)'3in- 


(3.S) 


gives  a  measure  of  the  error  for  the  maximum  shear  stfess  that 
is  associated  with  use  of  Eq .  (3.4).  The  ratio  given  by  Eq. 
(3.5)  is  plotted  in  Fig.  10  along  the  radial  line  9*20°. 

Equation  (3.2)  together  with  Eq.  (2.23)  can  be  rearranged 
to  the  following  more  convenient  form  for  the  purposes  of 
drawing  lines  of  constant  maximum  shear  stress  near  the  crack 


l. 

l^Tm^ 2 _  |~  (l-k)cos2a 

+  — T-  [(l-k^cos 
4(2)^  (L 


~i-1?||(l-k2)cos2a-(l-k)2cos22a]F1(9)+['ci-k)2sin4ajF2(e)|^ 

yj|[^(l  +  k)2-2(l-k2)cos2a+(l-k)2cos22ajF3(9)+  (1  -k)  2sin22ajF4  (0) 
+  ^4(l-k2)sin2ct-2(l-k)2sin4ajF5(0)|  -  0,  (3.6) 


where 


F^  (9)»sin6sin-y- 
F^  (6)*sin0cos4p+2sin|- 
F^ (9) *sin20 

For  a*7r/2,  Eqs.  (3.6)  reduce  to 


F4(0)=(4-3sin26) 

6 

Fs  (3)®sin0cos-j 


(3.7) 
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|sin26»0.  (3.8) 


The  pronounced  effect  that  the  load  biaxiality  has  on  the 
maximum  shear  stress  near  the  crack  tip  is  shown  in  Figures 
11-13.  It  is  observed  in  Fig.  11  that  as  the  load,  biaxiality 
passes  from  tension- tension  to  uniaxial  tension  to  tension- 
compression,  the  extent  of  high  level  shear  about  the  end  of 
the  crack  increases  appreciably.  Moreover,  since  plastic 
yield  is  strongly  influei  ced  by  the  maximum  shear  stress  (c.f. 
the  yield  conditions  of  Tresca  and  Mises,)  the  plastic  yield 
region  about  the  crack  tip  should  show  a  similar  pattern  of 
variation  with  load  biaxiality,  that  is,  should  increase  in 
size  as  load  biaxiality  varies  from  tension-tension  to  tension- 
compression. 


3.2  TWO  COLLINEAR  CRACKS 


Equation  (2.41)  combined  with  Eq.  (3.1)  produces  exactly 
the  same  form  for  the  square  of  the  maximum  shear  stress  as 
Eq.  (3.2)  for  the  single  crack,  except  that  and  ^  in  (3.2) 
is  now  replaced  by  K^j  and  Kjj  as  defined  by  Eq.  (2.42)  and 
(2.43).  When  these  expressions  for  the  stress  intensity 
factors  are  explicitly  introduced, the  equation  defining  the 
isostats  of  maximum  shear  stress,  analogous  to  Eq.  (3.6), 
becomes 


J 

Tm\  rClrk)co.s2a 
a  )  l  2 


fte) 


+_1T  l  .  (p )  j  [  (!•  k2 )  cos  2a-  ( 1  -k)  2cos 2  2a]  F 

\ij  ■  !s  ■ 


4(2) 


'h  Lj(P}j[ (l+K)2cos2a^(l-k)2cos22a] F3(0) 


.  2  ._2; 


+  [(i-k)"sin<-2a]F4(0)+((l-k‘*)sin2a-ii^lL  sin4a]Fs(0)  J  -'O',  j-1,2,  (3.9) 

where 

■  2 \  /•_ n  -2 

atb'-a")1 


bcp)  ■  ■ 


(3.10) 


(b  -a  ) 

and  the  F^(9),  k*l,2,...,5  are  as  defined  by  Eq.  (3.7).  Figures 
14-16  give  an  indication  of  the  variation  of  the  normalized 


•*$V. 


maximum  shear  stress  with  load  biaxiality  for  cracks  of  unit 
length  inclined  at  a=45°  and  spaced  at  2a=0.20. 


3.3  SHEAR  PANEL 

For  the  cracked  shear  panel  the  maximum  shear  stress  near 
the  crack  tip  is  determined  from  Eqs .  (2.56)  and  (3.1)  as 

Tm  3  8^:[K12sin29+K22(4-3sin20)+4K1K2sinecosyJ 

-  xsin2s| - — r  sin0sin-7p+ - — r-/sin0cos-^-+2sin4VTsin2 8 

LC2frr) 4  2  (2irr)*'  2  2/ 

(3.11) 

Again,  if  the  customary  singular  approximation  for  the  stress 
were  used  in  Eq . (3.1),  i.e.,  Eq.  (2.56)  without  the  2tsin23 
contribution,  the  result  would  be  only  the  first  line  of  Eq. 
(3.11).  Designating  the  first  line  of  Eq.  (3.11)  by  t*, 


-27- 


'm  8tt 


k[h 


2sin20+K,2(4»3sin20)+4K1K2sin6cos| 


]  ■ 


(3.12) 


a  measure  of  the  error  involved  by  the  use  of  Eq.  (3.12)  father 
than  Eq.  (3.11)  is  provided  by 


a-  32  j 


aiM  eos2  3(sin9sin-y-+2sin|-  j-sin2  3sin9sin4p 


[sin2  8sinz0+cot2  8cos2  3(4-3s.in49)-4cos2  8sin0cas 


-  8  I 


sin“ 8+001^2  8(4 -3s  in  0)  -4cot2  SsinQcos | 


(3.13) 


which  is  illustrated  by  Fig.  17. 

Equation  (3.11),  in  which  and  K2  are  given  by  Eq.  (2.58), 
can  be  written  as 

+  sin2B||)--j^  j  sin228F1(6)-sin2Bcos2SF2(0) 


-  | jsin228F3 (9)+cos 228F4 (8) -4sin28cos2 3F5 (9)  j  3  0  ,  (3.14) 

where  the  Fk(9)>  k*l , 2 , . . . , 5 ,  are  as  defined  by  Eq.  (3.7). 

Equation  (3  14)  may  be  used  to  draw  isostats  of  normalized, 
maximum  shear  stress  as  shown  in  Fig.  18, 
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SECTION  4 


ANGLE  OF  INITIAL  CRACK  EXTENSION 


The  maximum  tensile  stress  criterion  may  be  used  to  obtain 
an  estimate  of  the  angle  at  which  a  crack  extension  will  begin 
to  take  place.  Although  the  criterion  itself  cannot  be  considered 
to  be  wholly  satisfactory,  it  nevertheless  has  been  shown  to 
offer  rescnably  good  correlations  with  experimental  data  under 
tensile  loading  [9,10]  ,  and  can  serve  as  a  means  of  demonstrating 
the  influence  of  load  biaxiality  on  the  direction  of  crack 
extension. 

Designating  8Q  as  the  angle  which  locates  the  tangent  to 
the  direction  of  initial  crack  extension,  the  criterion  predicts 
crack  extension  along  the  radial  line  normal  to  the  direction 
of  the  maximum  tensile  stress  or,  equivalently,  along  the  dir¬ 
ection  parallel  to  the  plane  upon  which  the  tangential  stress 

component  t._  attains  maximum  value.  A  mathematical  statement 

y  w 

of  the  criterion  requires  that  at  some  arbitrarily  small  radial 
distance  0<r0<<l  from  the  crack  tip,  the  angle  of  crack  exten¬ 
sion  is  determined  by  simultaneous  satisfaction  of  the  con¬ 
ditions  : 


t09  (r0  ’  9(P  >0  ’ 


3tee(ro,90') 

Te~ 


8  t99  (r0,90^„ 

- T~  0 

30 


(4.1) 


The  unspecified  radial  dimension  rQ  is  necessitated  by  the 
singular  nature  of  the  stress  at  the  crack  tip  and  introduces, 
thereby,  a  small  measure  of  arbitrariness  in  the  predictions. 
It  should  also  be  observed  that  no  elastic  material  parameters 
appear  in  the  equations  for  stress  in  the  plane  theory  of 
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linear  elasticity  (neglecting  body  force),  consequently,  the 
above  criterion  predicts  the  same  angles  for  all  materials. 

A  series  expansion  for  the  tangential  stress  component 
t q q  can  be  shown  to  have  the  following  form  [3]  (up  to  terms 

r  ^ 

of  order  (-)  ) 

a 


s  A(r)  |B1(sin|+sin-y-)+B2(3cos|+cos^-)J +B3sin2e  (4.2) 

where  the  stress  t,  the  function  A(r)  and  the  coefficients 
B^,  B2,  B^,  are  defined  below  for  each  particular  load-crack 
geometry.  Insertion  of  Eq.(4.2)  into  the  conditions  (4.1) 
leads  to  the  set  of  requirements 


AfrjjjBj  sin-^.sin-Iij*B2^cos^cos-T^]J 
A(rg) |b2  3cos60-lj-3B^sine0J*4B 


<-cos-/Jj.B3sin2e0>  0 


3sin-j-cos6Q  *  0 


(4.3) 


0Or  2^0  e0  T  2®0*1) 

A(rQ)  B^sin-y  1  9sin^-j— 7  +  B2cos-y  6-9cos^-y-  j 


'  +  2B3cos26q  >  0 

which,  when  satisfied  collectively,  determine  the  crack  exten¬ 
sion  angle  9q. 


4.1  SINGLE  CRACK 


For  the  single  inclined  crack  the  undefined  quantities 
appearing  in  Eqs.  (4.2)  and  (4.3)  are  as  follows: 

t»a 


(k-  3)  sinacosot 
B2-k+(l-k)sin2a 
Bj* (l-k)cos2a 
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A  summary  of  the  variation  of  the  angle  of  crack  extension,  9Q, 
with  crack  inclination,  a,  and  load  biaxiality,  k,  is  shown 
by  Fig.  19,  for  which  the  arbitrary  value  of  Tq  was  chosen 
such  that  (a/r0)*0. 01 . 

4.2  TWO  COLLINEAR  CRACKS 

For  the  collinear  cracks,  the  undefined  quantities 

become 

t=*o 


A(r) "Aj  (r)  •  Lj  (p)  ,  j-1,2 


B^3 (k-  3) sinacosa 
B2*k+ (l-k)sin2a 
B3= (1 -k) cos2a 


Figures  20  and  <.i  illustrate  the  variation  of  crack  extension 
angle  with  load  biaxiality  and  inclination  of  the  cracks. 

Figure  20  is  for  t;  -  .nner  crack  tip  z1s*a  and  Fig.  21  for  the 
outer  crack  tip  22^0.  The  crack  spacing  was  taken  at  2aa0.20. 
According  to  Fig.  20,  the  inner  crack  ends,  although  closely 
spaced,  will  not  extend  along  their  original  directions  to  join 
together  when  the  cracks  are  inclined,  except  under  equal 


tension-tension  biaxial  load,  k*l,  or  when  the  cracks  are 
horizontal,  a*Tr/2,  and  k*2. 

4.3  SHEAR  PANEL 

For  the  cracked  shear  panel,  the  undefined  quantities 

become 


t  *  T 


B^*3cos2  0 
B2  =  sin2  8 


fc)'  0'01- 

The  variation  of  crack  extension  angle  with  crack  orientation 
is  shown  in  Fig.  22,  along  with  experimental  data  for  707S 
and  2024  aluminum  alloys  reported  by  Liu  [11]. 
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SECTION  S 

DISPLACEMENT  OF  THE  CRACK  BORDER,  COLLINEAR  CRACKS 


The  components  of  displacement  for  the  points  defining 
the  borders  of  either  one  of  the  inclined  collinear  cracks  can 
be  determined  using  Eq.  (2.36)  for  q(z)  and  w(z)  in  Eq.  ( 2 . 1)3 
for  the  points  along  the  x-axis  for  which  a<x<b.  Approaching 
the  crack  from  the  positive  side  one  obtains* 

2u[u<*>  (x,0)  +  iu^+)  (x,0)J-  (<-l)A^(  +  )  (x)-BX(+)  (x)J  -  (<+l)Cx, 


where 


4*  (x) 


f *  x2dx 
/  [(x2-a2)(x2-b2)] 


(5.1) 


jjbE(*,p)-|  (b2-x2)(x2-a2)jis 


(x)  =/X- 

i  1 0 


(5.2) 


x  “  -  a  " ) (x  -b  )  ] 


F  (d> ,p)  . 


F ( , p )  and  E ( d> , p )  are  the  elliptic  integrals  of  the  first  and 
second  kind,  respectively,  with  modulus  p.  Substitution  of 
Eqs .  (S.2)  and  (2.37)  into  (5.1)  gives  for  the  displacement 
components  along  the  positive  border  of  the  cut 


Approaching  the  cut  from  the  negative  side,  it  may  be  shown  that 
4^"^(x)  ■  -'F^(x),  X^'^(x)  3  -X*-  +  ^(x) 
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ux+)(x'0)“l:(1‘k)  j(<‘1)b  E^*p)i(|)2‘1][1'(l)2])!l  -*(?>(<!>, P)Jsin2a 

-  (k+1)  x  cos  2o }  '  /  a- 

2  1  (5‘3)  * 
Uy‘J(x,0)— ^  |(>-l)b  E(*,p)-g  ^  (|)  -l]  l-(f)  y  -»(p)F»,p)l[Cl»k) 

-  (l-k)cos2ot]-  (<+l)  (l-k)xsin2a  j  . 

For  horizontally  oriented  cracks  and 

x  x  y  y 

with  this  understanding 

Vx'0)  *  "k  (1'k)(K*1)  X  (5.4) 

•Uy(x.O)  •  -^;C<-l)bjEM,p)-f  ([(|)2-l][l-(|)2  • 

In  the  limit  a+0  leaving  a  single  crack  of  length  2b,  u  remains 

Jv 

the  same  and  Eq.  ( 5 . 4 ) 2  reduces  to 

uy(x,°)  *  (k-  1)  (b2*x2) **  ,  (5.5) 

which  coincides  with  the  result  obtained  from  the  single  crack 
solution.  Equation  (5.4)2  for  the  double  crack,  and  Eq.  (5.5) 

for  the  single  crack,  are  shown  in  Fig.  23,  each  plotted  for  a 

\  , 

crack  of  unit  length,  illustrating  the  effect  of  the  presence 
of  the  second  crack  on  the  uy  component  of  the  crack  border 
displacement  for  a  crack  spacing  2a-0.2.  The  influence  of  the 
horizontal  load  ka  on  the  crack  border  displacements  is  evident 
from  the  above  equations. 


) 
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SECTION  6 


GRIFFITH  CRACK  INSTABILITY  CRITERION 


The  criterion  for  fracture  introduced  by  Griffith  initially 
in  1921  [12],  and  restated  in  1924  [13],  is  expressed  as  follows: 

"The  fundamental  conception  of  the  new  theory 
is  this.  Just  as  in  a  liquid,  so  in  a  solid  the 
bounding  surfaces  possess  a  surface  tension  which 
implies  the  existence  of  a  corresponding  amount  of 
potential  energy.  If  owing  to  the  action  of  a  stress 
a  crack  is  formed,  or  a  pre-existing  crack  is  caused 
to  extend,  therefore,  a  quantity  of  energy  proportion¬ 
al  to  the  area  of  the  new  surface  must  be  added,  and 
the  condition  that  this  shall  be  possible  is  that  such 
addition  of  energy  shall  take  place  without  any  in¬ 
crease  in  the  total  potential  energy  of  the  system. 

This  means  that  the  increase  of  potential  energy  due 
to  the  surface  of  the  crack  must  be  balanced  by  the 
decrease  in  the  potential  of  the  strain  energy  and 
the  applied  forces." 

The  total  potential  energy  of  the  system  is  stipulated  to 
be  a  maximum  at  the  point  of  incipient  fracture  (unstable  crack 
extension) ,  or 


H  *  0  f6,1) 

at  the  onset  of  crack  instability.  When  applied  to  an  elastic 
ally  deformed  body  with  an  existing  crack  of  length  (a) ,  the 
total  potential  energy  of  the  system  consists  of  the  elastic 
potential  energy,  V,  and  the  crack  surface  energy,  r, 

P  -  V  +  r  ,  (6.2) 


where 


*/  ?tjkejkdV'(y*TkukdA  +  y"bkukdV)  *  )»kml»2»3  (6*3) 

R  'S  R  ' 


and 


U 


fk 


dV 


(6.4) 
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is  the  elastic  strain  energy  of  the  body,  which  is  subjected 
to  the  body  force  bk  and  the  surface  traction  applied  over 
the  region  R  bound  by  the  surface  S.  By  virtue  of  Clapeyron's 
theorem  [14]  however, 


f TkukdA  ♦ f  bkukdV  »  2U,  (6.5) 

S  R 


allowing  the  total  potential  energy  of  the  system  to  be 
equivalently  expressed  as 

P  -  T  -  U  ,  (6.6) 

with  the  corresponding  condition  for  crack  instability 

gj(r-U)  -  0.  (6.7) 

Since  Griffith  assumes  the  surface  energy  r  to  be  a  linear 
function  of  the  crack  surface  area  multiplied  by  a  constant 
surface  energy  density,  y,  application  of  the  fracture  cri¬ 
terion  is  reduced  to  determining  dU/da,  the  derivative  of  the 
elastic  strain  energy  of  the  entire  body  with  respect  to  the 
crack  dimension. 
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SECTION  7 


ELASTIC  STRAIN  ENERGY  DERIVATIVE 


An  estimation  of  the  load  that  must  be  applied  sufficient 
to  cause  an  existing  crack  in  a  structure  to  become  unstable, 
that  is,  to  suddenly  propagate  and  cause,  thereby,  complete 
separation  of  the  body,  can  be  obtained  from  application  of 
Griffith’s  fracture  criterion.  As  indicated  in  the  previous 
section,  all  that  is  necessary  in  order  to  apply  the  criterion 
is  knowledge  of  the  rate  at  which  the  elastic  strain  energy 
of  the  body  changes  with  extension  of  the  crack  size,  i.e.,  a 
calculation  of  the  derivative 


lim 

(a  '-a)-*-0 


(«)■ 


dU 

37 


(7.1) 

in  which  U  and  U'  represent  the  elastic  strain  energy  of  the 


body  having  the  crack  lengths  (2a)  and  (2a'),  respectively,  and 
where  the  crack  length  difference  Aa«a'-a  is  arbitrarily  small. 

In  the  case  of  two  collinear  cracks,  an  assumption  must 
be  made  at  the  outset  as  to  which  of  the  crack  tips  is  first 
anticipated  to  extend.  It  would  appear  on  the  basis  of  the 
ratio  of  the  crack-tip  stress  intensity  factors  at  the  inner 
and  outer  crack  ends,  Fig.  8,  that  for  relatively  closely 
spaced  cracks,  crack  growth  should  be  expected  first  at  the 
inner  crack  ends.  Thus,  corresponding  to  the  strain  energies 
(J  and  U'  are  the  initial  and  extended  crack  lengths  (b-a)  and 
(b-a ') ■b- (a- Aa) ,  respectively,  from  which 

/4TT 

(7.2) 


lim 

(b-a')-*- (b-a) 


U'-U 

TFTTn^iT 


-s 
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Consider  the  plane  body  with  a  single  horizontally 
oriented  crack,  (2a),  and  let  R  designate  an  arbitrary  region 
of  the  plane  bound  by  the  arbitrary  simple  closed  curves  C^, 

C2  that  enclose  the  crack  (of.  Fig.  24)..  For  the  plane  body 
with  horizontally  oriented .  col linear  cracks ,  (b-a) ,  the  region 
R  is  bound  by  three  arbitrary ^  s;iiiip l$y.ci6sed,, curves  *  k*l ,  2 , 3 , 

that  enclose  .  both,  cracks,  ■  (cf  Fig  .  2  5 ) .  yV-'; .  (• 

The  difference  of. elastic  strain  energy  (per  unit  thick- 
ness)  of  the  region;  R  of  the  body  having  incrementally /..differ¬ 
ing  crack  sizes  is  ;  :  ‘  v.  . 


U^CR)  -  U(R)  -JJ  i.(t  -kerk-tjkejk)dA,  ,  ,k-l  ,2.  (7.  3) 

From  the  linear  elastic  'constitutive  relations 

'  E  *t*jk 

the  following  identity  is  easily  established 

tjkejk  "  tjk®jk  *  ('tjk  +  tj'k)^j'k“ej;k)  “  ^tjkrtjk^ejk+ejk>  »  (7 

•  .  *  /  •  .  .  •  *  '  •  •  *  v  \ 

which  allows  Eq,  > (7 . 3 j  to.  be  given  the  equivalent  expressions 

' -  "  "  "  .  '  V  :: 

U"  (R)  -  U(R)  »  IjfjTvt  Jk.+^  j  k)  (ejk‘ej  k)dA  ,  j  Vk- 1,2,  '  (7,6) 


U"(R)-U(R)  *  l//  (tjk‘tjkMejk+ejk)ciA  >, 


(7.7) 


Either  expression  (7.6)  or  (7.7)  can  be  used  in  performing 
the  differentiation  indicated  in  Eqs .  (7.1)  and  (7.2). 


fc'jSi  j***'*’**  »*>••■  — 


Momentarily  writing  rj k+t j k** jk ’  ej  k'ej  k“ej*k  »  and 
uk’\*uk  ’  the  iirst.of  the  equivalent  expressions,  Eq.  (7.6), 


may  then  be  written  as 


■O'*  (R)  *  U(R)  * 


iff  £.*k  e*k aA  ■ 

R 

R 


(7.8) 


The  open  set  of  points  defining  the  domain  of  the  mathematical 
analysis^  which  consists  of  the /^tended  complex  plane  with  the 
points  Re (z) * j x | <a  defining  the  crack  removed  (or  for  the  two 
collinear  cracks  with  the  points' .|a.|<Re(z)*  |  x  |<  |  b  |  removed), 
is  a  multiply  connected  domain,  as  may  be  seen  by  recourse  to 
sterographic  projection  of  the  extended  cut  complex  plane  onto 
the  Riemann  sphere  [15 ] .  Thus  upon  application  Green's  theorem 
for  the  multiply  connected  domain  [16 j ,  taken  in  conjunction 
with  the  equations  of  equilibrium  and  the  assumption  of  zero 
body  force,  the  multiple  integral  given  in  Eq.  (7.8)  over  the 
region  R  may  be,  transformed  into  the  following  sum  of  line 

i  '  •  i  .• 

integrals  about  the  bounding  contours 


U-(R)-U(R) 


iE^[tLu**txyuy]dy-[tyyu/txyux]dx)  ’ 


(7.9) 


'J  '(R)  -U(R) 


iS^j[(^x+txx)(u;-ux)+(t;y+txy)(U;-uy)  ^ 


*[^t)T+t>y^lV'uy^+^txy+txy^ux'ux^jdx|  ’  C7*10) 


-39- 


where  N=  2  for  the  single  crack,,  and  N=3  for'  the  multiple  crack 
geometries .  ’V’:'.;  .•  ;  ‘V- : 

The  boundary  contours  are  arbitrary,,  hence  they  roa^ 


be  chosen  to 
may  be,  taken 


facilitate  computation.  The  outer  boundary , •  C.  ,  , ,, 


to  be  a  circle  centered  at  the  origin  with  an 
arbitrary  radius  r* • sufficiently.  large  so  as  .to  allow  the  - 
circle  to  completely  enclose  the  crack  (or  cracks).  The 
first  of  the  sum  of  line  integrals  in  EqV  (7,10)  then  has  the 
circular  line  integral  form 


:  <  i 


[u'.(R)-uCR)]c  -jTo'f\[ 

-  TT  '  . 


)(ux'ux) * < “y)j 


COS0 


r*rr 


m 

}tyy+tyy)  (‘uy"uy^+  ^txy+txy^  ^ 


(7.11) 


The  inner  boundary  (or  boundarie-s)  of  the  region  ,R  can  be  taken 
to  coincide  with  the  boundary  of  the  extended  crack  (or  cracks) . 


Accordingly,  the  second,  of  the  sum  of  line  integrals  in  ,£q <•  ;/,7 .,10) 

•  -  •=  ■'  ■  ■  •  ,s;'  >•'  ’ '  i’  ?y  ■■  ■ 

will  involve* integration  along  the  x-direction  only,  and' thus 


leave  the  reduced  forpt 


[u '  (R)  - U (P.)]  i ^  j  (t^y* tyy>  (u'-u 


y)*(^y+tiy)(^‘ux)|dx';  (7'12) 


‘2 


for  the  single' crack  ,  geometry .  and  '<%. 

,jV(RHJ(R) 


■  'V  i  'r  ' 
'  .  C2  C3 


(Sr  -*  <»•”) 


for  the  multiple  crack  geometry. 
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7.1  SINGLE  CRACK 

Further  evaluation  of  the  line  integrals  of  Eqs.  (7. 11) 
and  (7. 12)  requires  expressions  for  the  stress  and  displace¬ 
ment  i;.:  which  can  obtained froni'  £qs .  (2.11),  (2.12)  and  (2.1)*'. 
The  result ,  expressed  in  series  form,  is  as  follows: 


•xx ; 


j2V  (f):tco! 


[cos 26'  J  2sin6  sin30] 


/  \  4 

+  3 ( )  (j '  cos 40  -  s^n0  si'n50]  + 


yy 


7<2  +  {-)  (cos; 2 9  +  2sin0  sin30] 


, .  +  k(r)  { j  c o s 4 0  +  sine  sin59]  + 


(7.14) 


i  -i,  xy 


t  *  cr{(r).sin0  cos30  +  4-  (§1  sine  cos50  + 


ivr 


2  Vr 


and 


w' 


r  4  M  |2 

'  4 

+  j(  )[  3sin6sin49 


[<(l+k)  +  k-3]  r  cos0+  [sin0sin20 
4  v 


-(s1) 


cosfl] 


M 


a  (lr 

uy  3  4u|2f  K 


cos  30] +  . 

2 


(7. 15) 


(1+k) + (1- 3k) ] r  sin0+ 


ccs20] sin0 


4 . 

+  (j(<+l)sin39  -  3sin9  cos49]  + 

for  (a/r)<l. 

Along  the  circle  C1  with  r*r0,  it  follows  from  Eqs.  ( 7 . 1 4 T 


and  (7.15)  that 


2 

t'  +tvv  *  a  2k  +  ---{cos  2  6  -  2sin6  ,c  in39] 

'  .  .•*.  "  "  .  ro  ,  .  '  ■  ■ '  '  '  ■ '  ■ 

■f  K-  ^4  ~ j  cos46  -  Sin9  s inS el  + 


,  .  I'  '  '2  '  V 

tyy  +  tyy  *  a  2  +  |^-[cos29  +  2sin6  sin39]  ('7.16) 

1  ro  "  .  "...  .. 

'  +  2l*  ^  T  COS46  +  sine.  sinSe  ♦.  ....... 

'  •  r0  L  '  * 

Cxy+txy  3  sin0  c°s20  ♦  Sin9  cos5e  + 

I  r0  .  ■  ro 


Vux  *  ^pff£5i-n*  »i»M 


•(¥) 


cose] 


+  .  ~-ya(3sine  sin40  - 


cos  30 T .+ 


“j-'v  ■  ftjTfSs1)- 
•  ttRt) 


a3Aa 


cos2ej  sine 


sln36  -  3sin0  cos40 


(7.17) 


Note  that  in  the  above  expressions  terms  to  order  (Aa) ^  and 
righer  are  omitted.  This  entails  he  loss  of  accuracy  because 
Aa  may  be  taken  arbitrarily  small,  and  approaches  zero  in  the 
limit  process.  Substitution  of  Eqs .  (7.16)  and  (7.17)  into 
Eq .  (7.11)  leaves 
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(fi)-U(R) 


_  c/7  faf  I 


4aAa[k(sin0  sin20  cosQ 


■M 


COS^d) 


(7.18) 


(s1)-* 


2  2 
me-  cos26  sin  0] 


n=l x  u/ 


d0  . 


The  functions  H  (0;k),  nal,2,3,...,  appearing  in  the  integrand 
are  products  of  and  sums  of  sine  and  cosine  functions,  which 
are  finite  for  -ir<0<7r.  In  addition,  the  sequence  of  positive 
constants  { fa/ r q )  ^ }  approaches  zero  monotonical  ly  as  n-*-«>. 
Accordingly  for  each  n  and  -irsestr,  it  is  always  possible  to 
find  numbers  M  and  N  such  that  the  bound  functions  Hn(e;k) 

satisfy  |H^(0;k)+  .  +  Hn(0;k)|<M  for  any  n>N.  It  follows 

from  the  Dirichlet  test  for  convergence  [17]  that  the  infinite 
series  appearing  in  the  integrand  of  Eq .  (7.18)  converges 
uniformly  for  all  -ir*0*Tr.  The  order  of  integration  with 
summation  in  Eq .  (7.18)  may  therefore  be  interchanged,  allow¬ 
ing  it  to  be  reduced  to 


[u'.(R)-U(R)Jc  *  [k  (2  -k)  ♦  <+  2] 


n-l  u  -tt 


(7.19) 


For  the  infinite  sheet  geometry,  the  outer  boundary  contour 
Cj  of  R  must  approach  infinity.  In  the  limit  as  rQ-*-«>  the 
infinite  series  appearing  in  Eq .  (7.19)  converges  to  the  value 

.  .  /  *  \2n  /•■tf  _  .  I  /  _  \2  /  .  v4 


zero,  i.e. 


r“"-  £  fe)  Y>C6:kJde%*“'te)2H*tk,*(?ao)4  H2*tk) 


(7.20) 
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because  each  of  the  quantities  H^fk),  from  the  integrated 
bound  functions  H  (0;k),  is  a  finite  number  and  rn>a. 


0 


In  Eq.  (7.12)  C2represents  the  crack  border  with  crack  size 
2a  or  2a'*2(a  +  Aa).  Along  the  cut  |x|<|a|,  y*0,  the  stress  free 
surfaces  require  t  »t  »0.  Likewise  for  the  extended  crack  |x|<|a+Aa|, 

yy  xy 

y*0,  t'  -‘ mt'- =*0.  From  Eqs.  (7.14)  and  (7. IS)  it  is  clear  that  along  ya0 

yy  xy 


and  | x | < | a | ,  txy-uy»0. 


Because  of  the  symmetry  of  t  relative  to 


the  crack  plane  t^y^*t^y\  whereas  for  the  displacements  in  the  y- 


direction  u 


-(  +  )  -(-) 


The  positive  and  negative  superscripts 


y  y 

designate  the  opposite  faces  of  the  crack.  With  these  considerations 


in  mind,  Eq.  (7.12)  can  be  reduced  to 

'a+ Aa 


U  "  (R)  -  U  (R) 


JC2 


(7.21) 


By  means  of  Eqs.  (2.11),  (2.12)  and  (2.1), 


[  tyyJ  y  =  o 
x>a 


[uy]y-0 


-7  ? 

ax(x“-a^) 


x<a (a+Aa) 


^  [2aAa- (x2-a2)] 
4u 


so  that 


ju'CR)-O(R)]  c  -  -  j 

1  a 


a  +  AaUaAa-  (x2- a2) 

- r^i, — 

(  x  -a 


(7.22) 


(7.23) 


x  dx 


-  4a. 


(7.24) 


The  change  of  the  elastic  strain  energy  with  change  of 
crack  dimension  for  the  infinite  plane  body  can  be  obtained 
by  substitution  of  Eqs.  (7.19),  (7.20)  and  (7.23)  in  the 
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limits  (a  '-a.)  »Aa-*-0  and  rQ-*-°°: 

dU  _  lim  /u'-u\_  lim  [U' ( R) -U(R) ]  r, 
dl  (a'-a)^0VpT7/  Aa+O  Ai  k 

N  r0-»  k*l 

=  +  (7.25) 

The  derivation  of  Eq.  (7. 25)  proceeded  from  the  first  of 
the  two  equivalent  expressions,  Eqs .  (7.6)  and  (7.7),  for  the 
strain  energy  difference.  A  calculation  beginning  with  Eq. 
(7.7)  instead,  involves  only  a  difference  in  signs  for  the 
terms  appearing  in  the  integrands  of  the  line  integrals  of 
Eq.  (7.12).  In  place  of  terms  such  as  (t "  +t  )  (u'-u  )  there 

XX  XX  X  X 

now  appears  (t"  -t  ) (u/*+u  ) ,  and  so  on.  Stress  sums  are  re- 

XX  XX  X  X 

placed  by  stress  differences,  while  displacemnt  differences 
are  replaced  by  displacement  sums.  The  remaining  part  of  the 
calculation  is  identical  in  every  respect  to  the  one  above. 
The  result  is  the  same,  although  the  contributing  terms  are 
slightly  different,  i.e., 

d|  ,  (2-  k)  -<j  *  siojgZH  (K*1) 

.  L*V).  ^l*k(2-xr)j  .  (7.26) 

The  first  term  of  the  first  line  is  the  contribution  from  the 
outer  circular  boundary  in  the  limit  as  Tq^00  ,  while  the 
second  is  from  the  integral  about  the  boundary  C2  when  it 
coincides  with  the  crack  boundary. 
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If  one  can  imagine  crack  extension  during  which  the 
outer  boundary  of  the  body  remains  fixed  (a  so-called  'fixed- 
grips'  condition),  then  the  displacement  differences  (u'-u  ), 

XX 

(u'-u  )  are  zero  along  the  outer  boundary  and  the  contribution 

/  / 

to  the  strain  energy  difference  from  the  integration  about  the 
outer  boundary  vanishes.  This  leaves  only  t.he  second  term 
of  the  second  line  of  Eq .  (7.25)  for  the  strain  energy 
derivative,  which  comes  from  the  line  integration  about  the 
crack  border  boundary  C2 ,  i.e., 

*  ”  ^ (<+l)  <  0 ,  f°r  (<+l)  >  0, 

grips 

The  elastic  strain  energy  of  the  body  decreases  if  a  crack 
extends  under  'fixed-grips',  and  its  rate  of  decrease  is 
independent  of  the  load  biaxiality. 

If  one  can  imagine  crack  extension  during  which  the 
tractions  along  the  outer  boundary  of  the  body  remain  un¬ 
changed  (a  so-called  'dead-lead'  condition),  then  the  traction 
differences  Ct'x~txx) ,  (tyy‘tyy)  are  zero  along  the  outer 
boundary  and  the  contribution  to  the  strain  energy  difference 
from  the  integration  about  the  outer  boundary  vanishes. 

This  leaves  only  the  second  term  of  the  first  line  of  Eq.  (7.26) 
for  the  strain  energy  derivative,  which  also  comes  from  the 
line  integration  about  the  crack  border  boundary  C2>  i.e., 

(ai)dead  "  (,c+1)  >  °»  £or  (k*1)  >  °- 

load 
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The  elastic  strain  energy  of  the  body  increases  if  a  crack 
extends  under  ’dead-load',  and  its  rate  of  increase  is  in¬ 
dependent  of  the  load  biaxiality. 

These  two  quantities,  opposite  in  sign  but  identical  in 
magnitude,  are,  to  within  a  factor  of  two,  what  is  often 
identified  in  the  literature  as  Irwin's  "strain  energy  re¬ 
lease  rate"  G  [18] . 

The  right  side  of  Eq.  (7.2S)  expressing  the  strain  energy 
derivative  dU/da  can  be  viewed  as  a  continuous  function  of  the 
crack  size,  FCa;a,  k,  k) ,  with  a,  k  and  <  (or  v)  assuming  the 
role  of  parameters.  It  is  apparent  from  the  form  of  Eq.  (7.25) 
that  for  any  crack  size,  a,  and  applied  load,  a,  the  values 
assumed  by  the  strain  energy  derivative  may  be  positive  or 
negative  or  zero,  depending  on  the  values  of  the  parameters  k 
and  k.  The  latter  occurs  when 


(l+v)/(l-3v)  ,  for  plane  stress 


(7.27) 

l/(l-4v)  ,  for  plane  strain. 

Since  Griffith's  fracture  criterion  requires  solving  for  the 
extreme  value  of  the  potential  energy  by  setting  dP/da  ■ 
d(T-U)/da  *  0,  the  vanishing  of  one  of  the  two  contributions 
entering  into  this  expression  renders  the  calculation  meaning 
less.  Therefore,  for  the  condition  given  by  Eq.  (7,27) 
Griffith's  crack  instability  criterion  is  inapplicable. 
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7.2  TWO  COLLINEAR  CRACKS 

Because  of  the  symmetry  of  the  applied  loads  and  the 

•  \ 

crack  geometry  relative  to  the  x-axis  (horizontally  oriented 
cracks,  a*ir/2)  Eqs.  (2.1)  may  be  shown  [3,19]  to  simplify  to 

txx  a  Re(2*(z))  -  y  In  (2<T(z))  +  S 

tyy  «  Re(2$(z))  +  y  Im  (2$'(z))  -  S  (7.28) 

t xy  *  -y  Re  (2<J>"(z)) 

and 

2wu  *  (<-l)  Re(<Kz))  -  y  Im  (2<t(z))  +  S  x 

(7.29) 

2uuy  a  (k+1)  Im($(z))  -  y  Im  C  2  ( z ) )  -  S  y. 

The  constant  S  is  required  by  boundary  conditions  (2.6)  to 
have  the  value 

S  *  -  !  c(l-k)  .  (7.30) 

Equations  (2 . 35) - (2 . 37)  and  (7 . 28 ) - (7 . 30 )  lead  to  the 
following  series  expansions  for  the  stress  and  displacement 
components 

t  -  l^k  +  N (X)  Kr  (cos 2 0  -  2sin0  sin30) 

XX  L 


+  3M( A) 


cos49  -  sin9  sin59)+0 


71 


7(i 

t  *  y {"2  +  N ( A)  (cos29  ♦  2sin9  sin36) 

y>  *•  L 

+  3M  ( A )  — ^  1 cos40  +  sin0  s  in  5  9  ^  +  0  f 

r  '  'r  > 

txy  *  a[N(x)  ^7  sine  cos30  ♦  |M(A)  sine  cos50  +  °(^)j 


(7.31) 
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x  4  u 


U- 


(l  +  k) +  k -  3 )  r  cos8  + 


N(A)~-^sin6sin2  9-^^Jcosej 


+  j  M(X).^j(3sinesin48- |^yi)cos39)  +  0 


<$»] 


uy  *  &[?{ 


^||ic(l+k)+l-3k^r  sin0+N(X)^-^^ijsin 
j  M(X)^^^i)sin30-3sin6cos49^+O^jj 


in6'sin8cos20 


(7.32) 


where  r>b  and 


N(X)  -(fa/br  +  1  -  2X(p)] 


M  ( X)  »  £(a/b) 4  ♦  2/3(a/b)2  ♦  1  -  4/  3(^^.)x  (p)  .  (7.33) 

Since  the  crack  size  increment  (Aa)  is  arbitrarily  small, 
terms  of  order  (Aa)2  and  higher  may  be  ignored  in  writing 


cj  k  *  cjk  *  Hjk  (4a) 


uk  *  uk  *  ft*  (Aa)  - 


j  ,k«l ,2 


which  in  conjunction  with  Eqs.  (7.31)  and  (7.32)  lead  to  the 
following  expressions  for  the  integrand  of  the  line  integral 
in  Eq.  (7.11)  along  the  circle  C^,  r^rgJ 
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(txxnxx)“°j2k+[N(X)+lla  (Aa)]^T  (cos29*2sin6sin36) 

+  s|m(X)+j~  (Aa)j^ ^|cos48-sin6sin50^+O^.j j 

:)Vnyy)"0|2',‘[N(-X)+7  fa  (cos29+2sin6sin38) 

Jr0 

+  CAa)j^j^|cos40+sin6sin5e^+O^^J 

:xynxy)*°j[N(A)  +  fa  (Aa)K  sin6  cos3Q 

Jro 

+  j|m(A)  (Aa)J^-  sin8  cosSQ  +  ofer-M 


sin8  cos 58  +  01 


©l 


(7.3. 


x'ux)  *  (Aa)  7^  ^sine  sin29  '(iSr)c°se) 

+  J  If  C  Aa)^j(3sin8  sin46  '^ijccs3e)  ♦  °(^rjj 

y'uy  *  Tirjfl  (Aa)  7^({^r)sinQ  ‘  sine  cos29) 

+  t  fa  *•  Aa^  ^5‘^^y^^sin39-33in8cos40 


+  T  ff  C  Aa) 


in38-33in8cos48 


)*  '(S)l  ■ 


The  above  expression,  when  substituted  into  Eq.  (7.11), 


[u'(R)-U(R)]c  -  2^y  |k(sin6  cose  sin29  -(i^i)co528] 

“  7T  ' 

+[(^7^)sin29  *  sin40  cos2ej  ^  (Aa) 

+  Ua)2  (^-)2n  Tn  (9;F(p),E(p),k)  )  de  .  (7.36) 

n-1  0  / 

The  functions  T^ (0 ; F (p) ,e (p) , k)  include  products  of  sine  and 
cosine  functions  and  of  the  complete  elliptic  integrals.  As 
such  they  are  bound  function  for  0<a<b  and  -irsetir.  It  may  be 
shown  that  derivatives  of  N(X)  and  M(X) ,  which  involve  the 
elliptic  integrals,  again  produce  elliptic  integrals,  and  are, 
therefore, also  bound  Quantities  for  0<a<b.  Accordingly,  by 
the  same  argument  given  in  Section  (7.1),  the  series  converges 
uniformly  on  the  interval  [ -rr , tt]  and  it  follows  that 

[u'(R)-U(R)]Ci  ■  2^[k(2-K)  ♦  k  *  2]  H  CAS) 

(7.36a) 

*  2*ir  (4a)S  (h)j  V9;F(p)>Efp>k)  d«  • 

n-l  '  O'  * 

In  the  limit  rQ^«  the  infinite  series  appearing  in  the  above 
equation  converges  to  zero,  i.e., 

r0~  H  (rr)2n/'rTnce;F(p).H(p),k)de 

u  n-l  u 

(7.37) 

•ro-lft)7!^)^),”  +(^}4T2*(F(p)'E(P)‘k)+ . j"  0 
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be.c aus e  o'adh  of:  the  .quantities  Tn'  (F  (p)iE  (p) ,  k) ,  ..  from  the  in' 

V  te:grated  bound;;;funct.ions  Tn/(9;;  F  (p) ,  E  l.p)  ,:k)  •,  is  ;a  "finite  numbbr 
,.-p  I  n  the:  two  l.ine  •  .integrals  o  f  ,Eq  ;v  (  7  , 1'30.i.:  the  •  inner  "  - 

■  boundary  contours  C2  and  C,  coincide  with  the  borders  of  the 
c.r.ack.s For  the  extended  cracks5  the  boundary  conditions'  on 
-y  ! -  t he  crack  ..surf aces:,  requir^ythat, ' »  0  for'1' j  a>-ha  j<.  j:x  |,<  |  b  j.  \  . . 
whi  le  for  the  unextended  cracks  t  *  0  for  j  a  |<  |  x  j  <  |  b  | 

Because  of  symmetry  considerations  u^-0  for.  |  x  j<  |  a  1  ,  • 
uv  f+')  arid  ,  t^-t^  V  Consequently  ; '  Eq .  ..:(7r':3)  may  be  ' 


y  y  ■  yy  •  yy  ;-  .?s 

reduced  to ...; 


2  3  a " 


dt  .  (7 


.The  crack  -  tip  ■  expressions  Eqr. ,  C'2 . 41.)  2  and;  (2 . 47 )  ,  may  be 
employed  to  evaluate  Eq.  (7; 38)  since  the  interval  of  inte¬ 
gration  is  over  the  crack-tip  region.  Thus  for  the  stress 
just  ahead  of  the  unextended  crack-tip 


uyy] 


Ki(a> 

(  2  7T  t  ) 


h  ’ 


(7.39) 


while  for  the  displacement  just  beyond  the  extended  crack  tip 
[uvi  .  .  ’  — TJ - \%r)  t7-4°! 


r*Aa- t 


where  (a")*K1(a-Aa) .  Evaluation  of  Eq.  (7.38)  can  now 
proceed  to  give 


[u'(R)-U(R)"| 

L  j  2  3 


K:  (a) (a- Aa)  (<+l) 


(7.41) 
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Substitution  of  Eqs.  (7.36),  (7.37)  and  (7.41)  into  Eq.  ■; 
(7.2)  in  the  limits  as  (b-a'O+fb-a)  or  Aa-*0 ,  %  and  rg-*-rt,  yields 


dU;i  lim 

da  (b  -a ")  -*■  (b  -  a) 


aW 


(  U-  --u 
|  (F-T^'  -Tb'-aTT 


lim  ^  [lT(R)-U(R)lr, 

_  As  *■ 

Aa-*0 
rQ~° 


-  T1k(2!-K).K.21^ai  A-  Kx2  (a) 

‘ 1  •  I ,  -  • 

By  means  of  the  following  derivative  formulas 

■  ;  .  .  .  I . 

2" 


(7.4  2);.  . 


dF(p) 
da  - 


-  +  -^7  [E(p) 
ap 


F(P)] 


7 

dE(p)  *  +  [E (p)  -  F(p)J;  , 


3a 


(7.43) 


ap 


and  the  defining  Eqs.  (2.37)  and  (7.33)^,  it  follows  that 

A  2 


dN(A)  .  +  2IA(p)-(a/bp_ 
a [ 1  - (a/b) 2] 


<7; 44) 


From  Eq.  (2.43)^,  where  a=Tr/2,  and  Eq.  (7.44)  it  is  seen  that  , 


K2(a)  ■  o2*b2  l.HEj:.(a/bl,2I .  ,  .  Abi  *j(X)  . 

1  a  f 1 - (a/b) 2]  L  da 


(7. 45) 


Consequently, substitution  of  Eqs.  (7.44)  and  (7.45)  into  Eq 
(7.42)  leaves 


g  ,  ♦  g^Vk(2.K)1lMP?-(»/b^] 


4  u 


(7.46) 


a ( 1 -  (a/b)  ] 

It  will  be  recalled  that  calculation  of  Eq.  (7.46)  pro¬ 
ceeded  from  the  first  of  the  two  equivalent  expressions  (7.6) 
and  (7.7).  Were  the  derivation  to  have  followed  from  Eq.  (7.7) 
instead  of  Eq.  (7.6),  the  only  difference  would  involve  the 
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signs  appearing  in  the  integrands  of  the  line  integrals  in 
;  Hq.  (7. 10).,.  In  place  of  terms  such  as  (0+t/.„)  (u'-u  )  one 
wquld  have  (t .(¥x+ux^  arid  so  on,  i.e. .,  the  stress' sums 
.  are  replaced  t>y  stress,  differences  and  .  the, .displacement,-, 
differences  by  dispiac.ement  sums.  Ttie  remaining  part  of  the 
calculation  would  be  similar  'in  every  respect,  with,  the  result 


/.  ih  place  9 f  Eq .  (7 . 4  2fj . ;  'However  when  this  is:  combined  with 
( '’,♦44)  and  (-7. 45),  Eq.  (7.46)  is  reproduced. 

As  , til  the-  caS-e,  with:,  the  single  crack ,  the  derivative 
dU/da  given, by  Eq.  (7.46)  may  be  positive  or  negative  or'  zero, 
depending  on  the  values  of  the  parameters  k  and  The  latter 
possibility  -occurs  when,  the  load  biaxial i ty.  has  the  values 
given  by  Eq.  (7.27),  and  for  which  Griff ith 's ■ crack  instability 

■  i  i  " '  ,  '  .  ‘  •  '  . 

criterion  becomes  meaningless. 

' ■ '  ,  .  .<  •  .’  ' 

7 . 3  SHEAR  PANEL'  !  •  \  ' 

For  the  cracked  shear  panel,  (cf..  Fig.  9),  the  strain 
energy;  difference  expressed  by  Eq.  (7'.10)  for  the  bound  region 
R  shown  by  Fig.  26,  may  be  reduced  to 


tT(R)-U(R)»[U'(R)-U(R)Jr  +  (U-(R)-U(R)]  c 

L1  L2 


if.\ 

[<•; 


(t«*,o)tu;-V*(txy‘txy)(uy 


+  t  )  fu'-u  )  +  f t ■*  +t  )  fu 

yy  yy  y  y  xy  xy;  1 


co  s  9 


sin9  r^de 


/•a*4a  i(tf*>u'(t)»tc+)u'f>))dx 
J  2\  yy  y  xy  ux  / 


(7.4?) 


when  the  outer  boundary  C,  is  again  a  circle,  and  the  inner 
boundary  C£  coincides  with  the  borders  of  the  crack..'.'' 

The  stress  and  displacement,  obtained  t>om  Hqs . :  (2 . 50)>? 
(2.51)  and  (2.1),  are  given  by  : 

txx’l |2sin2 3+ (7)  t '2 (sin0cos30+sin20) cos2 3  ! 

+  2  (cos  20- 2sin0sin30)  sin2  3] +fx  (9  ;  0) 

.2 


(f) 


tyy53— j- 2sin2  3+^-J  [2sin0cos39cos28 


-  (2sin0sin39+cos2  8)  sin 2  8]  +fy  (0  ;  3)0^-^ 
txy*T|2cos2 3+ (r)  t (cos29"2sin0sin30)cos2S 


-2sin9cos30sin2  8]  +  fxy  (9;  8)0^ 


(7.43  ) 


and 
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•main  - 


UX  *  HI  j r (sin0cos23+cosesin2 3) 

:  i  a-  .  '  ■”  '  v-  v: .  -  'v  '' 

■!'*.  4  —  [  ( ( <+ 1 ) s in  8  +  2  s ine  c o  s  2 e } c qs 23  ;  ■  •-'/ 

V  +  { ( <  - 1 )'  CO s  9  •  2  s  ine s  in 2  0,}tsin  2 3]  +  ?x'f  9;;:3 )  0^2^  j  a:  ■ ' 

-  T  i  t.  (7,4.3) 

uy  3  217  {  r (cos  8cbs  2 3“.&'itt0S':in2 0-)*\ ^ C  { C k :> c o s B ;  '  V 
+2sin0sin29  f cos 23-  { (<+  i )  s in.e  -  2s in9cos  20 } sin 2.3 ]  '  ,  ' 


7-4 


■  i,;.' 


•  vv  Yw-:1  •'-''rp-v;"-;' :  .  ? :; 

' where-'  the  function  fx>  fy ,  gx  ,  and  gy  are  sine  and  cosine  , 
functions  of  the  ,variab2e  9,,and  th6  parameter  S,  , 

•  ...  From.  Eqs..  (7.4  8)  it:  may  be  seen  that  when  terms  of  higher 
'  or^ev  in  tIle  arbitrat'il.y  small  increment  of  crack  extension 

are  ignored,  then  along,  the  circular  boundary  Cn 

{  2  ;■  ./ 

( tyynyy)c  ~T)2sin2‘3  ;*  f ( 2.  $  inf  co  s  5  9  c  o  s  2  g 


J  •  •v,i'  +..(cos2e-2sine )  s;in2  3)  +fy  (6  ;  3) o^  j  j  ' 

^ *xx+txjP  **  j ^ c in23+--— 2  (si n6cos  3 0  -  s'in2 9 )  cos  2  8 

■'  '  <  ..  r0  ? 

-  (2&in9sin30+cos20)sin2S] +f^  (9;3)0|^4]j 

■'  Vr»;) 

*'txy  +  txy^ *T  |*cos 2S+^— (cos26  +  2sin0sin30)cos23 


+  2sin0cos30sin2  3]  +  f  (6  ;  S)0( 

'r0 


(7.50  ) 
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while  from  Eqs .  (7.49) 


(ux'ux) aJu  |  "r"~'  ^(<  +  l)  sin9-2sin0cos2e^  cos  23 

+  ^(<-l)cos0  +  2sin6sin20^  sin2  3j  +  gx(9; 

(Uy‘Uy)3j^  |  ~~  £^(K-l)cos9-2sin0sin29^cos2B 
(<  +  l)sin6  +  2sin9cos2  0\sin2  8l  +  gv(0;3)O 


(7.51) 


j  °y 

Substitution  of  Eqs.  (7.50)  and  (7.51)  into  the  first  of  the 
line  integrals  of  Eq.  (7.47)  yields 


2  C * 

[U"(R) -U(R) ]c  ~  J  (a • Aa (<  +  sin20-cos20 

^  -  7T 


+  (sin^2 3-cos42 3) (2cos29sin^9+sin^20) +sin43  (7.52) 

2  00 
+  sin4$  (2sin20sin40-3ssin49)  ]+a  •  Aa^ 

n»l 

The  functions  Gn(0;3)  are  sums  and  products  of  sine  and 
cosine  functions,  so  that,  as  before,  the  series  may  be 
shown  to  converge  uniformly  over  [  - rr , tt ] .  Thus 

Gn(9;S)de  ,  (7.53) 

where  it  is  clear  that  after  term  by  term  integration,  the 
resulting  infinite  sum  of  bound  quantities  will  converge  to 
zero  as  rQ-*-°°: 


U" (R) -U (R) 


2 

,  T  TTa 
4y 


2< • Aa+ 


2 

t  a 


“  /  „  \2n/-TT 

-Us)  /. 


ft) 


G  (0  ;  6)  }d0 


n 
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2n  r n 

lim  Zl  l  —  )  I  G  f  0 ;  3)d0  -  0  . 

V*  m  '  T0>  J  n  . 


(7.54) 


Because  the  crack  extension  increment  Aa  is  arbitrarily 
small,  the  second  of  the  integrals  of  Eq.  (7.47)  can  be 
evaluated  using  Eqs.  (2.56),  (2.57)  and  (2.58).  These  give 

■<);>]  e-o  -  -  X-§  ™  [‘S’]  ,-o  (7.SS) 


for  the  unextended  crack,  and 


-(  +  )' 

uy  0-7T 

‘  -* _ c 


t  /t  (a* Aa)  sin2  8(<+l )  JAa-x 
*  Zy~  1  2  tt 


r* (Aa-x) 


>Ct)]  9-w 


r= (Aa-x) 


+  Ju  ^<+1^ (x+a)cos2s 


t/tt  (a  +  Aa7cos2g(K+l)  JAa-x 

«-y-  \  2u 


+  TU  (x+a)sin2B 


(7.56) 


for  the  extended  crack.  The  second  integral  of  Eq.  (7.47) 
thus  gives 


IT(R)-U(R)  J 


'Aa-x 


(7.57) 


t  ira(<  +  l 


The  change  of  elastic  strain  energy  with  change  of 
crack  size  for  the  infinite  cracked  shear  panel  is  obtained 
from  Eqs.  (7.1  ),  (7.53),  (7.54)  and  (7.57)  in  the  limits 
rn-*-°°  and  Aa-*-0 ,  with  the  result  that 
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h 


dU  ,  lina  (2k)  -  (k+1)  -  t2tts 

da  4y  v  4y 


-j— —  (k-1).  (7.5S) 

It  is  noted  that  for  <*1,  or  for  v*l  in  plane  stress,  or 


for  v**i  in  plane  strain,  dU/da*0. 
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SECTION  8 


ALTERNATIVE  CALCULATION  OF  THE  ELASTIC 
STRAIN  ENERGY  DERIVATIVE 


A  determination  of  the  elastic  strain  energy  derivative, 
dU/da  ,  can  be  approached  in  a  somewhat  different  manner,  that 
is,  by  calculating  the  strain  energy,  U,  of  the  entire  cracked 
body  first,  and  then  obtaining  its  change  with  respect  to  the 
crack  dimension  by  a  formal  differentiation.. 

Referring  to  Fig.  (24)  for  the  single  crack,  Fig.  (25) 
for  the  collinear  cracks,  and  Fig.  (26)  for  the  shear  panel, 
the  elastic  strain  energy  (per  unit  thickness)  of  the  region 


R  is  specified  by 


*  7^jT[tXx13c'  +  tyy1y  *  txy 


(3u  3u  V 

~5y  +  ~l£/Jdxdy*  (8-1} 


Integration  over  the  plane  region  can  be  transformed  into 
line  integrations  about  the  boundary  curves  by  application 
of  Green's  theorem  as  discussed  in  the  previous  section: 


U(R) 


k*l  Lk 


+t  vuv]dy-[tvvu  vu  ldx)  ,  (8.2) 


xx  x  xy  y 


yy  y  xy  x 


where  N*2  for  the  single  crack  and  N»3  for  two  collinear 
cracks.  Again,  the  arbitrary  outer  boundary  C^  can  be  taken 
to  be  a  circle  of  arbitrarily  large  radius  rg,  while  the 
inner  boundary  (or  boundaries)  may  be  taken  to  coincide  with 
the  crack  borders.  For  the  single  crack  and  shear  panel 


geometries  such  choices  will  give 
U(R) 


(8.3) 


(8.4) 


'l  '  T  ro/  lltxxY,xy'V'„coseMtyyVVx1  sine!de 

*  TT  T^0  r*r0  ' 


h  ’ 

and 


[tyyuy  +  txyux^dx  *  0  * 


(8.5) 


U(R)  -  lx  *  I2  ♦  r3 


Ix  *  Eq.  (8.4) 

I,  *  -  \{<P  *  [tyyuy  *  txyux]dx  -  0 


(8.6) 


(8.8) 


for  the  collinear  crack  geometry.  The  integrals  I,  about 
the  crack  borders  vanish  because  of  the  stress  free  crack 
boundary  conditions. 


8.1  SINGLE  CRACK 

Utilizing  Eqs.  (7.14)  and  (7. IS)  along  the  circle  rarQ, 
Eq.  (8.4)  gives 

2  00 

Il"TOr/|rOF(0{k)+a2GC0ik)*a2Z  (^)2n  Hn(6;k)  j  de  ,  (8.9) 

in  which 

F(e;k)»k(<(l+k)+k-3]cos29+[<(l+k)+(l-3k)]sin“9 
G ( 9  ;  k)*2k[sin9cos9sin29--y  (tc-l)cos29  ] 

+  (<(l+k)+(l-3k)]cos39cos9sin29+[(<+l)-2cos26]sin29 
+  y( < (l  +  k)  +  (k- 3) ] [cos29-2sinesin39 1  cos  2 9 

+  y[<(l+k)+(l-3k) ] [cos29+2sin0sin39]sin“9 
+  [<(l+k)+ (k-3) ]cos9cos39sin29  , 
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which  are  products  and 


and  the  functions  H  (0;k),  n*l,2,3... 
sums  of  sine  and  cosine  functions,  are  bound  over  the  interval 
-TTsesTr.  The  infinite  series  in  the  integrand  of  Eq.  (8.9) 
is  similar  in  form  to  the  series  appearing  in  the  integrand 
of  Eq.  (7.18),  and  may  be  proved  to  converge  uniformly  over 
the  interval  of  integration  by  the  same  argument.  Consequently, 
integration  of  the  first  two  terms  of  the  integrand  followed 
by  interchange  of  the  order  of  integration  with  summation  of 
the  series  leaves 


The  infinite  series  in  Eq.  (8.10)  vanishes  in  the  limit  rp-*-00 
in  the  same  manner  as  shown  by  Eq.  (7.20). 

By  means  of  Eqs,  (8.3)  -  (8. 5)  and  Eq.  (8.10),  the  elastic 
strain  energy  for  the  entire  infinite  body  is  obtained  in  the 
limit  rp-*-<»: 

U  »  lim  (i.)  +  l7  , 

Tp-*™  \  1/  2 

with  the  result 
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The  first  term,  UQ,  which  is  independent  of  the  crack  dimen¬ 
sion,  represents  the  strain  energy  of  the  infinite  non-cracked 
body  with  boundary  tractions  applied  as  in  Fig.  1.  The  second 
term,  which  includes  the  crack  dimension,  indicates  that  the 
elastic  strain  energy  of  an  infinite  sheet  of  unit  thickness 
containing  a  centrally  located  line  crack  is  different  from 
that  of  a  continuous  sheet  by  the  amount  AU.  It  is,  therefore, 
often  referred  to  in  the  literature  as  the  energy  of  the 
crack.  Formal  differentiation  of  Eqs.  (8.11)  and  (8.12)  with 
respect  to  the  crack  dimension  gives 

af  *  a?  (VAU>  ’  -’IF*1  i  .  (8.13) 

which  concides  with  Eq.  (7. 25). 

It  should  be  noted  that  for  equal  tension- tension  loads, 
k*l,  uU  given  by  Eq,  (8.12)  is  identical  to  Griffith's  result 
of  1921  [12],  and  to  Swedlow's  calculation  in  1965  [20].  For 
uniaxial  load,  k*0,  Eq.  (8.12)  reduces  to  the  expression 
obtained  by  Wolf  [21]. 

8.2  TWO  COLLINEAR  CRACKS 

Substitution  of  Eqs.  (7.31)  and  (7.32)  into  (8.7)  produces 
2  C  ^  ”  2n 

IimfeilroF(0;k)+N(A)b2G(e»k)+b2L(f)  Hn(9;A(p)’k)  |de  *  (s.14) 

^  n*l  G 

in  which 
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F (6 ; k)  -  kf<(l+k)+(k-3) ]cos0+[<(l+k)+(l-3k) ]sin  9 


G(9;k)  *  2k(sin9sin29  -  cose  ]  cos9  +  [  <  (l+k)  +  (1  -  3k)  ]cos9cos39sin20 

+  j[<(l+k)+(k-3)] [cos29-2sin0sin38]cos20+[ (<+l) - 2cos26 ] sin2 6 
+  i[<(l+kj+(l-3k) ][cos20+2sin0sin30]sin29 
+  [<(l+k)  +  (k-3) ]cos9cos39sin29 

The  functions  H  (6;X(p)fk)  are  products  and  sums  of  sine  and 
cosine  lunctions  and  the  complete  elliptic  integrals,  and 
are  bound  functions  over  -7rc0sir  and  0<a<b.  The  infinite  series 


in  Eq.  (8.14)  is  similar  in  form  to  the  series  appearing  in 
Eq.  (7.36),  and  may  therefore  be  proved  to  converge  uniformly 
over  the  interval  of  integration  by  the  same  argument.  Con¬ 
sequently,  integration  of  the  first  two  terms  of  the  integrand 
followed  by  interchange  of  the  order  of  integration  with 


summation  of  the  series  leaves 


2  2 

.  a  *r0 


o  2  irb  2 


[(<+l)k  +2(<-3)k+<+l]  +  -  ~  N ( X) [ l+k  (2  - <) 1 


X (p) , k)  de 


(8.15) 


The  infinite  series  in  Eq.  (8.15)  vanishes  in  the  limit  as 
r0-*-°°  in  the  same  manner  as  shown  by  Eq.  (7.37). 

By  means  of  Eqs.  (8.6)-(8.8),  Eq.  (8.14)  and  Eq.  (8.15), 
the  elastic  strain  energy  for  the  infinite  body  follows  in 
the  limit  as  rg^*: 

U  *  rii!  (!l)  *  !2  *  *3 


giving 


where 


U  -  UQ  4  AU  , 


(8.16) 


lim  l*2*r0 


0  V" 


YSf  [  (<+l)k  +  2 (<- 3) k  *  k  *  1] 


md 


(8.17) 


2  2 

AU  -  N (X)  1 1  ♦  k (2 * <)  ] 


In  these  expression  Uq  and  AU  have  the  same  interpretation 
as  given  in  Section  8.1  for  the  single  crack.  In  the  limit 
as  a-*-0,  leaving  a  single  crack  of  length  2b,  N(X)-*-l  and  AU 
assumes  the  value  appropriate  to  a  single  crack,  i.e.,  Eq. 

( 8  . 1 2 ) 2 •  For  crack  lengths  (b-a)»£  held  fixed  while  the 
crack  spacing  2a-*-<»,  leaving  the  infinite  body  without  cracks, 
N ( X) -*-0  and  U»UQ.  A  formal  differentiation  of  Eqs.  (8.16)  and 
(8.17),  utilizing  Eq.  (7.44),  reproduces  Eq .  (7.46) 

~  *  zz  cvau)  *  +  4^(i+k(2-<))i^l^^41-  •  c8-18) 

da  da  0  4M  a (1-  (a/b) 2  ] 


8.3  SHEAR  PANEL 

Substitution  of  Eqs.  (7.43)  and  (7.49),  with  rTg,  into 
Eq.  (8.4)  gives 


-IT 


r2(cos22B4sin22B) 


4  F  ( 0  ;  3 )  +  a 


2E(^)2"  Gn<9^]  dS 

n-1  u  J 


(8.19) 


) 


-65- 


where 

F(6;3)=(l-4cos223)sin220- (4sin29cos40 
2 

+2cos  0+l)cos20-2 (2sin8sin30+s in49sin4 3) +< 
and  G  (0;S),  n*l,2,...,  represent  functions  composed  of 
products  and  powers  of  sine  and  cosine  functions.  Integra¬ 
tion  of  the  first  two  terms  of  Eq .  (8. 19). leaves 
r  „  t2na2  ,  -2-2 '  ,  -  x2n 

h  Zu  *  —fir  (< 


-•  v. 


8)d0  .  (8.20) 


-ir  n*l 

The  infinite  sum  in  Eq.  (8 . 20,),  as  .  ijj  Eq.  (7.15),  converges 
uniformly  over  --rrsesir  by  the  same  argument  given  previously. 
Thus  upon  interchange  of  integration  and  summation,  it 
follows  that  in  the  limit  r^-*00 


r0-<* 


a  \2nf"  ) 

* Z 

■t)  j 

<  n*l  X 

0  7  <.  _  J 

7T  7 

I  a2[(tf  G*(6)  *(t)4  G2(B) 


(8.21) 


since  each  of  the  numbers  Gn(3)  are  bound  quantities. 

From  Eqs.  (8.5)-(8.5),  Eq .  (8.20)  and  Eq.  (8.21),  the 
elastic  strain  energy  for  the  infinite  cracked  shear  panel 
is  determined  in  the  limit  as  -r^^i 


v-  ('.) 


from  which 


U  -  a0  +  AU 


(8.22) 
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where 


u  ,  iikP*"T0  V 

■  vo^A^r/ 


arid 


AU 


T2  2 
t  Tra 

“§T“ 


(8.23) 


(<-l) 


..  A’.. 


Formal  differentiation  of  Eqs.  (8.22)  and  (8.23)  recovers 
Eq.  (7.56)  '  ,•••  ' ''  .  . 


dU 

..3a 


ZUU0  + 


2„  2 

r  ?ra 

4y 


(<-l) 


(8.24) 
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SECTION  9 


FRACTURE  LOAD 

The  tensile  load  necessary  to  cause  onset  of  crack  in¬ 
stability  (fast  crack  propagation)  according  to  the  Griffith 
crack. instability  criterion  may  be  determined  from  Eq .  (6.7j. 
The  surface  energy  of  the  crack  (per  unit  thickness)  is 
assumed  to  be  equal  to  a  constant  surface  energy  density,  y, 
multiplied  by  the  surface  area  of  the  crack. 

T  *  4ya  .  (9.1) 


9.1  SINGLE  CRACK 

For  the  single  (horizontally  oriented)  crack,  Eqs.  (7.25), 
(9.1)  and  (6.7)  lead  to  the  critical  applied  tensile  stress 
values 


1 8Ey 
I  rra 


(3v-l)k+ (l+v) 


and 


8F.Y 

r  i  ij 

*a(l+v!) 

[_(4v- 1 )k+l  J 

plane  stress 


plane  strain. 


(9.2) 


The  variation  of  the  fracture  load  with  load  biaxiality  and 
Poisson’s  ratio  for  plane  stress  and  plane  strain  is  shown 
by  Figs.  27  and  28.  Referring  to  Fig.  27,  for  Poisson  ratio 
values  v<l/3,  ac  increases  with  increasing  load  biaxiality, 
while  for  values  greater  than  one-third  it  decreases  as  the 
load  increases.  In  the  range  of  values  0.3<v<0.4,  which  en¬ 
compasses  most  structural  materials,  the  biaxial  effect  on 
the  fracture  load,  while  present,  is  rather  small  for  the  k 
variation  -lsksl. 


-68- 


For  values  of  k  equal  to  those  given  by  Eq.  (7.27),  the 
unbound  value  for  a indicated  by  Eq.  (9.2)  should  not  be 
interpreted  as  meaning  that  these  equations  predict  an 
infinite  fracture  stress.  As  remarked  in  Section  7,  for 
the  plane  infinite  centrally-cracked  body  at  these  particular 
values  of  the  load  biaxiality,  dU/da  is  zero  and  Griffith's 
fracture  criterion  becomes  inapplicable. 


9.2  TWO  COLL INEAR  CRACKS 


For  the  (horizontally  oriented)  collinear  cracks,  Eqs. 
(7.46),  (9.1)  and  (6.7)  give 


*  {a[l-  (a/bTl }  16yy 


& 


X(p)-Ca/b) 


r 


irb 


r  rr^rr-TU 


(9.3) 


The  variation  of  with  load  biaxiality  k  in  the  case  of 
plane  stress  is  shown  by  Fig.  29.  For  materials  with  Poisson’s 
ratio  less  than  one-third,  the  breaking  stress  must  be  in¬ 
creased  as  tensile  loads  of  increasing  magnitude  are  applied 
parallel  to  the  cracks.  This  trend  is  reversed  when  the 
Poisson's  ratio  value  is  greater  than  one-third.  The  curves 
shown  in  Fig.  29  for  the  multiple-crack  are  very  similar  to 
those  determined  for  the  single  crack,  i.e.,  Fig.  27.  The 
comments  given  above  concerning  load  biaxiality  values  at 
k*l/ (<-2) for  the  single  crack,  apply  in  the  same  way  for 
the  collinear  cracks. 


To  gain  some  idea  of  the  influence  of  the  crack  spacing 
on  the  breaking  load,  b  can  be  assumed  to  be  fixed  with  the 
value  b*l,  so  that  the  pair  of  cracks  have  the  fixed  lengths 


£=(l-a).  The  corresponding  critical  normal  tensile  load  for 
a  single  crack  of  length  (1-a)  is 


The  ratio 


!c  .  (9.5) 

°C  (2)^[ X  (p) -a^  ] 

provides  a  measure  of  the  degree  to  which  the  body  is  weak¬ 
ened  by  the  presence  of  the  second  crack.  The  variation  of 
Eq .  (9.5)  with  the  ratio  of  the  distance  between  the  cracks 
to  the  crack  length,  i.e.,  p=(2a/Jl),  is  shown  in  Fig.  30, 
from  which  the  following  important  observation  can  be  made: 

As  the  spacing  between  the  col  linear  cracks  approaches  four 
times  the  individual  crack  length,  the  critical  load  foi-  the 
doubly-cracked  body  is  just  about  the  same  as  the  critical 
load  of  the  body  with  a  single  crack  of  the  same  length. 

Thus  a  pair  of  widely  spaced  collinear  cracks  in  a  relatively 
large  body  behave,  with  respect  to  instability,  as  if  the 
body  contained  only  a  single  crack  having  the  length  of  one 
of  the  collinear  cracks. 

The  critical  load  for  collinear  cracks,  each  of  length 
£*(l-a),  can  be  compared  to  the  critical  load  for  a  single 
crack  having  the  greater  length,  £=2>2(l-a),  by  the  ratio 


(ad-a2)^ 

A(p)-a2 


(9.6) 
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c.f.  Fig,  31.  For  cracks  relatively  closely  spaced,  the 
critical  normal  tensile  load  necessary  to  cause  the  inner 
crack  ends  to  extend  and  join  exceeds  the  critical  load  for 
the  resulting  longer  crack.  Fracture  of  the  collinear  cracks 
will  therefore  proceed  in  two  stages:  The  inner  crack  ends 
will  extend  first  to  join  up  forming  the  larger  single  crack, 
whereupon  the  sustained  load  level  will  be  sufficient  to 
cause  the  second  stage,  onset  of  fracture  from  the  ends  of 
the  newly  formed  single  crack.  From  these  observations  the 
assumption  of  crack  instability  for  the  doubly-cracked  body 
in  terms  of  initial  extension  of  the  inner  crack  ends  is 
j  ustified. 

9.3  SHEAR  PANEL 

The  critical  applied  tensile  stress  for  the  cracked 
shear  panel  has,  by  virtue  of  Eqs,  (7.53),  (9.1)  and  (6.7), 
the  value 

8EY  ,q  . 

ma (1+v) (k- 1)  *  (  •  ) 

which  is  seen  to  be  independent  of  the  angle  of  crack 

orientation.  Note  that  for  k- 1 ,  or  equivalently,  for  v*l 

in  plane  stress  or  in  plane  strain,  dU/da*0  and  the 

above  expression  is  meaningless,  that  is,  Griffith's  crack 

instability  hypothosis  is  no  longer  applicable. 
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SECTION  10 


EFFECT  OF  LOAD  BIAXIALITY  ON 
FATIGUE  CRACK  GROWTH  RATES 

The  phenomenon  of  crack  formation  and  the  subsequent 
incremental  extension  of  such  cracks  under  repeated  or  cyclic 
load  application  until  catastrophic  fracture  ensues,  may  be 
usefully  thought  of  as  three  hypothetically  distinct  or 
separate  physical  events:  (i)  micro-crack  nucleation,  ter¬ 
minated  by  the  appearance  of  macro-cracks,  followed  by  (ii) 
crack  'propagation'  (which  is  here  taken  to  mean  very  minute 
incremental  crack  extensions  caused  by  cyclic  load  application)  , 
and  finally  at  the  critical  crack  size  (iii)  onset  of  rapid 
crack  propagation  to  fracture.  This  somewhat  artificial 
division  has  proven  to  be  fruitful,  nevertheless,  because  it 
allows  f r r  the  introduction  of  fracture  mechanics  concepts 
that  have  been  able  to  correlate  a  wide  range  of  fatigue 
crack  growth  experimental  data  appropriate  to  category  (ii) 
semi-empirically,  in  a  more  or  less  satisfactory  manner. 

10.1  SEMI -EMPIRICAL  FATIGUE  CRACK  PROPAGATION  LAWS 

According  to  the  generally  accepted  micro  structural 
theory  concerning  fatigue  crack  development  and  growth,  the 
basic  mechanism  of  crack  nucleation  and  crack  extension  under 
cyclic  load  is  cyclic  slip  (micro-plasticity),  with  slip  bands 
(multiple  slip  planes).  Crystalline  slip  results  from  the 
motions  of  dislocation  which  are  generated  by  the  action  of 
shear  stresses  along  potential  slip  planes.  Fatigue  micro 
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crack  nucleation,  which  leads  to  the  development  of  a 
macroscopic  crack  by  coalesence  of  many  mi.cro-cracks ,  and 
the  subsequent  propagation  of  this  crack  under  continued 
cyclic  load  application,  appears,  therefore,  to  be  intimately 
related  to  the  extent  of  plastic  deformation  that  developes 
around  the  crack  tip. 

This  realization  has  served  as  the  physical  basis  for 
the  development  of  a  variety  of  semi - empirical  fatigue 
crack  propagation  laws,  evolving  ultimately  into  the  widely 
used  form  [22-26] 

^  *  B (AKX ) n  ,  (10.1) 


where  (a)  is  the  current  crack  length,  N  the  number  of  load 
cycles  and  AK1*  CK1)rnax“  a  measure  of  the  applied  load 

range.  On  a  log-log  graph  this  relation  plots  as  a  straight 
line 


3M  *  4,118  +  n 


(10.2) 


for  which  B  represents  the  AK^«1  intercept  and  n  the  slope 

of  the  straight  line.  Relation  (10.1)  has  been  shown  to 

correlate  a  broad  range  of  fatigue  crack  growth  data  in  the 

moderate  to  high  aK^  range,  where  the  data  was  obtained  mostly 

from  tests  performed  on  center  cracked  sheets  (with  horizontal 

crack  orientation)  cycled  in  uniaxial  tension  between  (K.) 

l'max 

and  CK1)min*0 ,  that  is,  for  the  load  ratio  ^l^min/  (K^) 


max 


3min/a  v  *0 . 
max 
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Over  the  full  range  of  crack  growth  rates,  however, 

-  3  -  9 

ranging  from  10  to  10  inches  per  load  cycle,  fatigue 
crack  growth  rates  are  not  observed  to  fall  along  a  straight 
line  (cf.  Fig.  32).  At  the  higher  load  levels,  as  the  fatigue 
crack  approaches  the  critical  crack  size,  the  increment  of 
crack  growth  per  cycle  of  load  increases  causing  an  upward 
inflection  away  from  the  straight  line,  and  characterizes 
the  passage  from  the  fatigue  crack  propagation  stage  of 
category  (ii)  to  the  fast  fracture  stage  of  category  (iii). 
Here  as  the  curve  asymptotically  approaches  the  upper  vertical 
line  (da/dN)-*«.  At  the  lower  end,  the  downward  inflection 
from  the  straight  line  represents  progressively  smaller 
crack  growth  increments  per  lead  cycle,  so  small  as  to  be 
considered  negligible,  i.e.,  (da/dhO^O,  The  asymptotic 


approach  to  the  lower  vertical  line  is  therefore  interpreted 
as  defining  a  fatigue  threshold  value,  AK^ ,  below  which 
fatigue  crack  propagation  cannot  take  place  [27], 

To  take  account  of  the  frequently  observed  effect  of 
the  load  ratio,  E,  on  the  growth  rates,  [28],  as  well  as  to 
incorporate  the  departure  from  linearity  of  the  upper  end  of 
the  growth  rate  curve,  Eq.  (10.1)  may  be  modified  in  the 


manner  proposed  by  Forman  [29] 


C(Ak\  )m 

da  ,  _ l _ 

ZN  (l-UjK  -akx 


(10.3) 


where  the  parameters  C  and  m  are  determined  in  the  same  way 
as  the  parameters  B  and  n  of  Eq.  (10.1).  As  (Ki)max^Kc» 
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indicating  the  onset  of  fast  fracture,  the  denominator  of 
(10.3)  approaches  2ero  and  (da/dN)-*-».  The  deviation  from 
linearity  at  both  ends  of  the  spectrum  of  fatigue  crack 
growth  rate  can  be  accounted  for  by  the  form  proposed  by 
Hartman  and  Schijve  [301. 


da 

3N  ' 


m 


(10.4) 


(1-R)KC-AK1 

According  to  this  expression,  as  the  value  of  the  load  range 
AK^  approaches  the  lower  limit  or  fatigue  threshold  value, 
AKq,  (da/dN)-*0.  This  is  represented  in  Fig.  32  by  the 
asymptotic  approach  to  the  vertical  at  negligibly  small 
crack  growth  rates  approximating  (da/dN^O. 


10.2  EFFECT  OF  LOAD  BIAXILITY 

Any  survey  of  the  large  body  of  existing  fatigue  crack 
growth  data  identifies  many  factors  that  can  possibly  effect 
fatigue  crack  propagation  [28,  31].  These  may  be  grouped 
into  the  following  broad  categories: 

(i)  Geometry  -  including  crack  size,  shape  and  orientation, 

specimen  thickness  and  relative  size. 

(ii)  Loading  -  including  load  magnitude,  max-min.  load  ratio, 

biaxiality,  sequence,  cyclic  frequency. 

(iii)  Material  Properties  -  including  the  elastic  modulus, 

uniaxial  yield  and  ultimate  stresses,  cyclic 
stress-strain  properties. 

(iv)  Environment. 
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Such  a  large  number  of  possibilities  probably  explains 
why  for  any  given  one  of  them, one  almost  invariably  encounters 
conflicting  experimental  data.  Some  of  the  factors  no  doubt 
must  mitigate,  or  entirely  nullify,  their  effects  on  fatigue 
crack  growth  rate  when  considered  jointly.  Without  the 
benefit  of  a  comprehensive  rational  theory  of  fatigue  crack 
propagation  that  can  explicitly  account  for  many  of  these 
factors  and, thereby,  provide  a  basis  for  weighing  their  relative 
contributing  effects,  development  of  a  full  understanding  of 
fatigue  crack  phenomena  must  continually  entail  partial 
understanding  and  accompanying  confusion. 

With  respect  to  the  influence  that  load  biaxiality  might 
have  on  fatigue  crack  propagation,  qualitatively  important 
insights  and  expectations  may  be  gotten  from  the  information 
provided  in  the  previous  sections.  The  belief  that  the  elas¬ 
tic  stress  intensity  factor  (assuming  here  a  horizontal 
crack  orientation)  is  sufficient  to  completely  determine  the 
stress  about  the  crack-tip,  is  shown  by  equations  (2.24) 
and  (3.3)  to  be  incorrect  in  situations  where  loads  are  bi- 
axially  applied.  Figure  11  demonstrates  the  pronounced 
effect  that  the  horizontal  load  ak  has  on  the  intensity  and 
spatial  distribution  of  the  maximum  shear  stress  which,  in 
turn,  must  cause  a  qualitatively  similar  pattern  of  dependence 
on  ak  of  the  plastic  deformation  around  the  tip  of  the  crack. 
Thus,  the  pattern  of  variation  of  the  lines  of  maximum  shear 
stress  with  increasing  tension- tension  load  biaxiality  implies 
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correspondingly  reduced  plastic  deformation  and,  therefore, 
a  correspondingly  reduced  fatigue  crack  growth  rate. 

Equation  (9.2)^  and  Fig.  27  indicate  that  for  materials 
with  a  Poisson  ratio  in  the  vicinity  of  0.30,  the  fracture 
stress  ac  increases  approximately  linearly  with  increase  of 
tensile  load  biaxiality  in  the  load  range  - 1 < k<  2.5.  The 
plane  stress  fracture  toughness  for  a  horizontally  oriented 
center- cracked  sheet  of  finite  width  should  therefore  also 
show  a  similar  dependency  on  k,  i.e., 

Kc(k)=ac(k)  (7Ta)lsf(a/w)-[ac(0)  +  3k]  (ira) 4f (a/w)  (10.5) 

where  3  is  the  slope  of  the  a vs.  k  curve  of  Fig.  27  and 

f(a/w)  is  the  correction  factor  for  the  finite  width,  w,  of 

the  test  specimen.  A  similar  argument  can  be  made  for  the 

plane  strain  fracture  toughness  K  based  on  Eq.  (9.2) -  and 

1C  1 

Fig.  23.  The  semi - empirical  expressions  (10. 3) - (10. 5)  thus 
appear  to  indicate  that  an  increase  in  the  tensile  load 
biaxiality  should  cause  a  reduction  of  the  fatigue  crack 
growth  rate  relative  to  the  growth  rate  for  uniaxially  cycled 
load . 

The  reduced  crack-tip  plasticity  associated  with 
increasing  tension-tension  load  biaxiality  might  possibly 
diminish  the  effect  of  crack  closure  that  is  believed  to 
accompany  crack  propagation  through  cycically  loaded  and 
unloaded  plastically  deformed  regions  [27],  Any  such  miti¬ 
gation  or  lessening  of  their  effect  would  have  the  consequence 
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of  increasing  the  crack  growth  rate,  offsetting  thereby,  to 
some  extent,  the  anticipated  crack  growth  rate  decrease  with 
horizontal  tensile  load  increase  indicated  above. 
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PART  II.  EXPERIMENT 


The  agenda  for  the  experimental  part  of  the  research 
program  was  determined  by  selection  of  those  tests  that  could 
provide  directly  measureable,  or  observeable,  information  on 
the  influence  of  load  biaxiality  on  conditions  near  the 
crack-tip  and  on  fracture.  The  tests  that  were  performed 
involved  crack-tip  photoelasticity,  measurements  of  angles  of 
initial  crack  extension,  fracture  load  and  fatigue  crack 
propagation  rates,  all  under  varying  levels  of  horizontal 
load  application. 
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SECTION  11 


DEVELOPMENT  OF  THE  BIAXIAL  TEST  FACILITY 


To  be  able  to  perform  the  indicated  tests,  it  was 
necessary  at  the  outset  to  design  and  assemble  a  biaxial 
test  specimen,  and  a  testing  system  compatible  with  it. 
Several  of  the  more  important  features  of  this  part  of  the 
test  program  are  described  initially. 

11.1  SPECIMEN  DESIGN 

The  preliminary  biaxial  test  specimen  design  was  based 
upon  a  specimen  configuration  developed  by  Radon  and  co¬ 
workers,  [32],  which  was  shown  to  provide  a  relatively  large 
uniform  stress  field  for  the  given  specimen  size.  Several 
other  designs,  [33-37],  were  also  examined,  which  suggested 
several  minor  changes  in  the  preliminary  design.  A  photo¬ 
elastic  study  was  employed  to  ensure  that  the  design  would 
be  adequate  to  provide  the  desired  uniform  stress  field 
through  out  a  large  central  portion  of  the  center  cracked 
specimen  over  a  wide  range  of  the  biaxial  load  ratio.  The 
specimen  configuration  that  was  finally  adopted  is  shown 
in  Fig.  33.  A  second  large  specimen  was  also  designed  in 
the  event  that  a  larger  working  area  became  necessary.  This 
alternate  design  is  shown  in  Fig.  34.  It  was  also  noted  from 
this  initial  photoelastic  study  that  the  biaxiality  of  the 
applied  load  has  a  clearly  observeable  effect  on  the  crack- 
tip  stress  field  that  is  in  qualitative  agreement  with  the 
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predictions  of  the  analysis.  These  results  will  be  pre¬ 
sented  and  discussed  in  a  following  section. 

It  was  also  necessary  to  design  and  purchase  (or 
fabricate)  grips  to  interface  between  the  test  system  and 
the  specimen.  The  specially  designed  grips  that  were  pur¬ 
chased  performed  satisfactorily  throughout  the  testing 
program. 

11.2  TEST  SYSTEM  DEVELOPMENT 

The  second  phase  involved  designing  a  biaxial  test 
system,  which  was  comprised  of  an  existing  MTS  testing  system 
of  100  kips  capacity  as  the  vertical  axis,  a  portable  test 
frame  of  50  kip  capacity  as  the  horizontal  axis,  and  a  control 
system  capable  of  controlling  both  axis.  Since  only  two 
independent  actuators  were  planned  for  this  test  system 
(rather  than  four),  the  center  of  the  test  specimen  would 
not  remain  stationary  but  would  move  in  response  to  the 
motion  of  each  actuator.  In  order  to  prevent  this  motion 
from  introducing  side  loads  on  the  specimen,  the  horizontal 
axis  was  suspended  in  its  working  position  by  elastic  ropes. 

In  this  manner  the  test  frame  could  move  relatively  freely 
with  the  center  of  the  specimen.  However,  the  motion  of  the 
frame  under  rapid  fatigue  loading  conditions  could  introduce 
dynamic  loads  which  would  be  transmitted  to  the  specimen. 

The  magnitudes  of  these  loads  were  calculated  for  the  ex¬ 
pected  worst  case  conditions  and  were  found  to  be  small 
enough  to  be  neglected. 
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The  tost  system  design  incorporated  considerable 
versatility  so  that  a  wide  range  of  tests  could  be  per¬ 
formed,  The  control  system  purchased  as  part  of  the  second 
channel  (horizontal  axis)  was  designed  and  connected  so  that 
it  could  control  both  actuators  for  applying  phased  static 
or  cyclic  loads,  with  the  phase  angle  readily  adjustable  to 
any  desired  value.  Since  the  specimen  could  serve  as  a  load 
path  between  the  two  control  channels,  a  specially  designed 
pair  of  service  manifolds  was  incorporated  into  the  syst  m, 
effectively  eliminating  "cross-talk"  between  the  channels. 
These  service  manifolds  also  permit  completely  independent 
operation  of  the  two  channels  so  that  different  types  of  load 
programs  could  be  applied  to  each  channel.  For  example,  a 
randomly  varying  cyclic  load  could  be  applied  to  one  axis 
while  a  constant  amplitude  fatigue  load  or  a  static  load 
could  be  applied  to  the  other.  The  two  systems  could  also 
be  operated  independently  when  biaxial  tests  were  not  being 
performed. 

In  order  that  the  horizontal  axis  be  lightweight  and 
relatively  easy  to  manage,  the  horizontal  load  frame  was 
designed  and  manufactured  from  aluminum  structural  sections 
and  high-strength  aluminum  plates.  The  load  frame  design 
was  completed  and  the  components  were  fabricated  in  the 
school  machine  shop.  The  assembled  test  system  weighed 
less  than  1000  pounds,  of  which  the  actuator  alone  weighed 
approximately  600  pounds.  After  assembly  of  the  horizontal 
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test  system,  this  system  was  operated  for  a  brief  check-out  perio 
and  then  interconnected  with  the  existing  test  system.  The 
combined  system  was  then  checked  out  and  found  to  operate 
as  designed.  A  schematic  drawing  of  the  combined  test 
system  is  seen  in  Fig.  35,  and  photographs  of  the  completed 
system  are  shown  in  Fig.  36. 
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SECTION  12 


EFFECT  OF  BIAXIAL  LOADS  ON  THE  CRACK-TIP 
ISOCHROMATIC  PATTERNS 

As  the  biaxial  test  specimen  design  and  the  biaxial 
test  system  was  being  developed,  an  initial  photoelastic 
study  was  performed  for  the  purpose  of:  (1)  examining 
photoelastically  several  aspects  of  the  test  pecimen  design, 
and  (2)  providing  preliminary  information  about  the  possible 
effects  of  biaxial  load  on  the  crack-tip  stress  field. 

The  photoelastic  study  sought  to  answer  several  ques¬ 
tions  that  concerned  important  details  needed  to  complete 
the  design  of  the  test  specimen.  The  first  parameter 
examined  was  the  radius  of  the  fillet  between  two  adjacent 
load  tangs.  Initial  biaxial  specimen  designs  developed 
by  other  researchers  employed  large  circular  radii  or  hyper¬ 
bolic  curves  between  the  load  tangs.  Consequently ,  a  large 
radius  was  employed  initially  for  the  first  photoelastic 
specimens.  In  this  series  of  experiments  it  was  seen  that 
the  biaxial  forces  did  have  a  significant  effect  on  the  crack- 
tip  isochromatic  pattern,  however,  the  quantitative  agreement 
between  the  theoretically  predicted  isochromatic  patterns 
and  the  photoelastic  results  was  not  very  good.  Subsequent 
study  showed  that  the  large  fillet  radius  introduced  inter¬ 
actions  between  the  stresses  parallel  and  perpendicular  to 
the  crack.  Other  studies  (38]  have  confirmed  this  condition, 
and  have  also  shown  that  there  was  significant  nonuniformity 
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of  the  stress  field  along  the  width  of  the  tang.  Additional 
studies  by  other  researchers  have  been  made  into  the  effects 
of  different  geometrical  details  at  the  intersection  of  two 
adjacent  tangs  [39] .  The  effects  of  different,  small  radii, 
and  of  undercutting  were  examined  and  it  was  concluded  that 
employing  a  small  fillet,  as  in  Fig.  33,  provided  the 
optimum  configuration  in  the  intersection  region. 

The  second  parameter  examined  was  the  crack  length- to- 
specimen  width  (a/w)  ratio.  The  biaxial  photoelastic  speci¬ 
mens  were  fabricated  initially  with  small  crack  lengths  and 
after  photoelastic  examination,  the  cracks  were  extended 
and  subsequent  photoelastic  examinations  performed.  The 
results  of  these  studies  showed  that  the  extent  of  the  biaxial 
effect  was  not  altered  by  changes  in  the  (a/w)  ratio.  As  ex¬ 
pected,  the  crack-tip  stress  fields  were  affected  by  the 
(a/w)  ratio,  with  the  fringe  patterns  increasing  in  size  as 
the  crack  length  was  increased.  However,  the  effect  of  load 
biaxiality  on  the  fringe  patterns  was  not  significantly 
affected  by  changes  in  the  (a/w)  ratio. 

The  third  parameter  to  be  studied  was  the  crack-tip 
radius,  with  two  different  notch  radii,  and  a  natural  crack 
employed  in  the  study.  The  notch  radii  were  filed  using 
screw-head  files  of  different . radii ,  while  the  natural 
crack  was  produced  by  means  of  a  xazor  blade.  As  with  the 
(a/w)  study,  it  was  determined  that  the  sharpness  of  the  crack- 
tip  radius  had  a  significant  effect  on  the  crack-tip  stress 
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fields.  However,  the  magnitude  of  the  load  biaxial  it/ 
effect  was  not  measurably  influenced  by  variation  of  the 
crack- tip  radius. 

Since  the  effects  of  several  important  test  parameters 
were  shown  to  have  insignificant  interaction  with  the  pre¬ 
sence  of  the  horizontal  loads,  a  study  of  the  biaxial  load 
effects  on  the  crack-tip  maximum  shear  stress  was  then  per¬ 
formed.  The  most  extensive  data  were  taken  from  one  of  the 
early  photoelastic  specimens  having  an  (a/wj  ratio  of  approxi¬ 
mately  0.5,  a  natural  crack-tip,  and  a  large  fillet  radius. 

The  variation  of  the  lines  of  maximum  shear  stress  with  load 
biaxiality  indicated  by  the  analysis  of  Section  3.1,  and 
summarized  in  Fig.  11,  provides  a  basis  for  comparison  be¬ 
tween  the  theoretically  predicted  patterns  and  the  photo- 
elastically  observed  patterns.  The  observed  patterns  are 
shown  in  Fig.  37  for  k=0,  0.5  and  1.0.  A  comparison  of 
Figs.  11  and  37  show  good  agreement  for  the  angular  shift  in 
the  stress  profiles,  and  careful  measurements  confirmed 
the  agreement.  The  change  in  the  size  of  the  isostatic 
loops  with  increasing  tensile  load  biaxiality  agreed  with 
the  predicted  pattern,  although  the  magnitudes  did  not 
initially  correlate  very  well.  However,  when  a  correction 
factor  appropriate  to  the  large  fillet  radius  was  applied, 
the  sizes  of  the  isostatic  loops  and  the  associated  stress 
intensity  factors  agreed  much  more  clrsely  with  the  theoretical 
predictions . 
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SECTION  13 

FRACTURE  TESTING  OF  SPECIMENS  WITH  HORIZONTAL  CRACKS 

Upon  completion  of  the  test  system  development,  test 
specimen  design  and  photoelastic  verification,  the  fracture 
testing  program  was  initiated.  It  was  determined  that  three 
initial  test  series  would  be  performed  and  approximately  forty 
specimens  were  prepared  for  testing.  Two  of  the  test  series 
were  performed  on  707S-T6  aluminum,  0.063  inch  thick,  with 
the  center  slot  oriented  perpendicular  to  the  rolling  direc* 
tion  for  one  series  and  parallel  to  the  rolling  direction  for 
the  other.  The  material  chosen  for  the  third  test  series 
was  0.25  inch  thick  plexiglass,  because  its  material  pro¬ 
perties  were  quite  different  from  7075-T6. 

Thin  center  slots  were  machined  in  all  specimens  and  the 
slot  tips  were  subsequently  sharpened  to  a  tip  radius  of 
0.003-0.005  inch  by  use  of  a  screw  head  file.  The  7075-TG 
specimens  were  tested  with  the  slot  tips  in  this  condition 
(no  fatigue  precracking)  since  prioT  uniaxial  testing  ex¬ 
perience  with  this  material  indicated  that  subcritical  crack 
growth  occurs  prior  to  unstable  fracture.  With  the  existence 
of  subcritical  crack  growth,  unstable  fractures  initiated 
from  the  existing  crack  and,  therefore,  were  not  influenced 
by  the  notch  tip  radius.  In  contrast,  the  plexiglass  speci¬ 
mens  fractured  in  a  manner  apparently  more  brittle  than  the 
7075-T6  specimens.  It  was  anticipated  that  variations  in 
the  crack  tip  radius  could  cause  significant  differences  in 
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the  fracture  behavior  of  the  plexiglass.  Consequently, 
natural  cracks  were  introduced  into  these  specimens  by  press* 
ing  the  slot  ends  with  a  suitably  reinforced  razor  blade. 

The  final  crack  length-to-specimen  width  ratio  was  between 
0.4  and  0.5  for  all  of  these  tests. 

Since  707S-T6  is  a  relatively  brittle  aluminum  alloy  it 
was  decided  to  perform  two  additional  test  series  on  2024-T3 
aluminum.  The  interest  in  performing  the  2024-T3  tests  was 
based  on  two  considerations.  First,  2024-T3  is  an  important 
material  for  aerospace  applications,  and  it  is  often  treated 
as  a  prototype  material  for  the  2xxx  series  of  aluminum  alloys. 
Second,  the  additional  ductility  provided  by  2024-T3  was 
desirable  because  finite  element  calculations  by  others 
suggest  that  the  biaxial  effect  on  fracture  load  is  increased 
with  increased  material  plasticity.  These  specimens  were  made 
according  to  the  same  design  and  were  also  machined  from 
0.063  inch  sheets.  A  total  of  fifteen  tests  were  performed 
on  2024-T3,  in  which  the  starter  slot  was  parallel  to  the 
roll  direction  for  seven  tests  and  perpendicular  to  the  roll 
direction  for  eight  tests. 

In  performing  these  fracture  tests,  ramp  loads  were 
applied  to  the  two  perpendicular  pairs  of  specimen  tangs  in 
porportion  to  the  desired  biaxiality  ratio,  k.  For  the  7075- 
T6  specimens  the  biaxiality  ratios  varied  in  the  range 
0<k<1.8  and  for  the  plexiglass  specimens  the  range  was 
-0.5sk*2.0.  For  the  tests  on  7075-T6  and  2024-T3  the  loading 
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was  interrupted  periodically  for  the  purpose  of  recording 
the  loads  and  displacements  on  both  axis  with  the  use  of  a 
datalogger.  In  addition,  the  loads  and  displacements  per¬ 
pendicular  to  the  slot  were  recorded  continuously  on  an  x-y 
plotter.  For  the  plexiglass  tests  the  loading  was  not 
interrupted,  due  to  the  time-rate  effects  associated  with 
this  material.  Only  the  loads  and  displacements  perpendicular 
to  the  crack  were  recorded  on  the  plotter, 

13.1  FRACTURE  TEST  RESULTS  FOR  707S-T6 

For  the  two  707S-T6  test  series,  the  effect  of  load 
biaxiality  on  the  fracture  load  (perpendicular  to  the  plane 
of  the  starter  slot)  is  seen  in  Figs.  38  and  39.  The  speci¬ 
men  thickness,  B-0.063  in.,  was  considered  to  be  sufficiently 
small  to  establish  a  plane  stress  situation.  This  assump¬ 
tion  was  confirmed  by  the  mode  of  fracture,  in  which  the 
fracture  surfaces  were  oriented  at  45°  to  the  lateral  sur- 

i 

face  of  the  specimen  (100  percent  oblique  fracture) .  For 
both  test  series  the  critical  fracture  load  increased  15-20 
percent  as  k  was  increased  from  0  to  1.8.  Higher  k  values 
could  not  be  tested  successfully,  because  the  horizontal  tabs 
pulled  off  before  unstable  crack  propagation,  thus  releasing 
the  biaxial  constraint.  Since  the  crack  length  and  specimen 
dimensions  were  identical  for  all  specimens,  the  critical 
stress  and  Kc  values  increase  proportionally.  The  linear 
regression  fit  to  the  data  also  showed  the  same  increasing 
trend  in  both  test  series,  with  the  slope  of  the  lines 
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exhibiting  the  same  general  uptrend.  Poisson's  ratio  for  this 
material  is  approximately  0.30.  For  this  value  of  Poisson's 
ratio  Fig.  27  predicts  a  small  increase  of  the  fracture  load, 
which  is  confirmed  by  the  experimental  results.  The  difference 
in  results  between  Figs.  38  and  39  is  attributed  to  ehe  aniso¬ 
tropy  introduced  by  the  rolling  process. 

13.2  FRACTURE  TEST  RESULTS  FOR  PLEXIGLASS 

The  predicted  variation  of  critical  fracture  load  (or 
Klc)  with  k  shown  in  Fig.  28  suggests  that  failure  under 
plane  strain  conditions  (square  fracture  surfaces),  combined 
with  values  of  Poisson's  ratio  greater  than  0.2S,  should 
produce  the  reverse  trend  of  decreasing  critical  load  (or 
rlcJ  with  increasing  k.  Therefore  a  test  series  was  per¬ 
formed  on  plexiglass  specimens  with  B«0.2S  in.,  a  thickness 
which  was  large  enough  to  cause  plane  strain  fracture. 

Poisson's  ratio  for  this  material  is  high,  0.40<v<0.47, 
so  that  in  view  of  the  test  data  trend  obtained  for  the  7075 
aluminum  alloy,  a  definite  decreasing  trend  in  the  Klc  data 
would  provide  a  strong  confirmation  of  the  results  of  the 
theoretical  analysis.  The  results  of  the  plexiglass  tests 
are  shown  in  Fig.  40,  and  do  exhibit  a  trend  of  decreasing 
fracture  toughness  with  increasing  load  biaxiality.  These 
results  are  presented  in  terms  of  critical  Kj  rather  than 
critical  load  since  the  length  of  the  natural  cracks  was 
not  the  same  for  all  specimens.  Although  there  was  a  greater 
amount  of  scatter  in  this  data,  the  linear  regression  confirms 
the  decreasing  trend. 


Additional  confirmation  of  the  plexiglass  data  can  be 
found  in  the  test  results  recently  published  by  Radon,  et  al, 
(40 ,  41],  although  they  did  not  come  to  the  same  conclusion, 
possibly  due  to  the  lack  of  a  linear  regression  analysis  of 
their  data.  The  results  of  forty  tests  on  polymethylmeth¬ 
acrylate  (PMMA) ,  which  is  essentially  the  same  material  as 
plexiglass,  are  shown  in  Fig.  41.  The  best  fit  straight  line 
was  added  to  their  reported  data,  and  it  is  seen  that  the  two 
test  series  agree  remarkably  well  with  regard  to  the  magni¬ 
tudes,  the  trends,  (decreasing  slopes),  and  the  scatter  in 
the  toughness  values.  Although  the  extent  of  test  data  re¬ 
ported  here  is  insufficient  for  definitive  conclusions, 
nevertheless,  the  qualitative  agreement  with  the  predicted 
trends  for  both  materials  offers  considerable  evidence  as 
to  the  presence  of  biaxial  load  effects  on  the  fracture 
stresses  and  the  associated  K  values. 

13.3  FRACTURE  TEST  RESULTS  FOR  2024 -T3 

In  order  to  correlate  as  closely  as  possible  with  the 
7075-T6  tests,  one  series  of  seven  specimens  was  prepared 
with  the  starter  slot  parallel  to  the  rolling  direction  and 
a  second  series  of  eight  specimens  was  prepared  with  the 
starter  slot  perpendicular  to  the  rolling  direction.  The 
specimen  geometry,  the  slot  length  and  notch  tip  radius  and 
the  test  procedure  were  the  same  as  for  the  7075-T6  tests. 

The  effect  of  the  biaxiality  ratio,  k,  on  the  fracture 
load  for  2024-T3  is  shown  in  Figs.  42  and  43  for  the  starter 
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slot  oriented  parallel  and  perpendicular  respectively  to  the 
rolling  direction  of  the  material.  The  behavior  of  these 
two  test  series  is  clearly  self  consistent  in  that  they  both 
display  a  peak  in  the  critical  fracture  load  corresponding 
to  k  values  around  0.5  to  0.7S.  The  rate  of  variation  of 
critical  load  with  k  is  greater  in  the  range  of  0<k<0.75  than 
was  evidenced  by  the  7075-T6  results,  and  suggests  that  the 
increased  biaxial  effect  may  be  due  to  the  greater  ductility 
of  the  material.  However,  the  reasons  for  the  existence  of 
a  peak  value  and  subsequent  decrease  in  c  ritical  load  with 
increasing  k  is  not  apparent.  One  possible  explanation  for 
this  behavior  is  that  the  increased  ductility  of  the  material 
was  responsible  for  the  tangs  parallel  to  the  starter  slot 
to  fail  first,  thus  eliminating  the  biaxial  constraint. 

This  behavior  would  lead  to  the  prediction  that,  at  the  higher 
k  values,  the  breaking  load  would  approach  the  result  for 
uniaxial  tension,  k=0.  It  is  seen  that  both  tests  series  do 
display  this  type  of  behavior,  with  the  exception  of  the  test 
at  k*1.25  in  Fig.  42.  There  is  some  further  experimental 
evidence  to  support  this  type  of  behavior,  although  these 
observations  need  considerably  more  study  before  reliable 
conclusions  can  be  drawn.  The  difference  in  critical  load 
between  the  two  test  series  is  due  to  the  anisotropy  effect 
and  is  of  the  same  order  of  magnitude  as  the  effect  for  7075- 
T6 . 


SECTION  14 


FRACTURE  OF  ANGLE-CRACKED  SPECIMENS  OF  7075-T6  ALUMINUM 

Since  a  crack  in  a  structure  subjected  to  biaxial  applied 
loads  will  not  normally  be  aligned  with  either  of  the  princi¬ 
pal  loading  directions,  knowledge  of  biaxial  load  effects  on 
angle-cracked  specimens  is  of  considerable  practical  impor¬ 
tance.  The  loads  on  both  axis  of  the  angle  cracked  specimens 
must  be  considered  because  the  loads  on  each  axis  can  be  re¬ 
solved  into  components  parallel  and  perpendicular  to  the 
starter  slot.  The  component  perpendicular  to  the  starter  slot 
represents  the  principal  force  causing  fracture,  while  the 
other  component  coincides  with  the  horizontal  (biaxial)  load  for 
the  flat  crack  geometiy. 

Two  test  series  were  performed  on  7075-T6  aluminum  sheets, 

0.063  in.  thick,  with  the  starter  slot  oriented  45  degrees  to 
either  of  the  load  axis.  Four  tests  were  performed  in  each 
test  series,  with  one  test  at  each  of  the  biax^ality  ratios  of 
k  -  0,  0.25,  0.5  and  1.0.  In  one  series  the  primary  (larger 
load)  axis  was  aligned  with  the  rolling  direction,  while  for  the  othe 
the  primary  axis  was  perpendicular  to  the  rolling  direction. 

Figure  44  shows  the  broken  specimens  from  the  equal  biaxial 
tests  (k»1.0)  for  both  test  series,  where  the  principal  load 
direction  is  indicated  by  the  arrow. 


The  specimen  and  starter  slot  were  prepared  in  the  same 
manner  as  for  the  specimens  having  the  starter  slots  oriented 
perpendicular  to  the  primary  load  axis.  The  same  testing 
procedure  employed  previously  for  the  aluminum  test  series 
was  also  followed  for  these  tests.  The  load  versus  displace¬ 
ment  for  the  principal  axis  was  plotted  continuously  on  an 
x-y  recorder,  and  the  loading  ramp  was  interrupted  periodi¬ 
cally  for  the  purpose  of  recording  the  load  and  displacement 
values  for  each  axis  using  a  datalogger.  The  critical  or 
maximum  load  in  the  principal  loading  direction  was  also 
recorded  using  a  digital  multimeter  with  memory  capability. 

According  to  the  analysis  of  the  angle-cracked  specimen 
given  in  Section  2,  the  opening-mode  stress  intensity  factor 
is  given  by  [cf.  Eq.  (2.23)] 
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( (1  +  k) * (1 -k)cos2a ] , 


corresponding  to  an  infinite  sheet  geometry,  where  the  angle 
of  crack  inclination  a  is  defined  in  Fig,  1.  For  these  tests 
a=450.  In  addition,  a  finite-width  correction  factor  must  be 
appended  to  the  above  expression, which  then  becomes 

K  *  ^~Al+k)f(^)  .  (14.1) 

If  the  critical  value  of  the  applied  stress,  obtained  as 
described  above,  is  inserted  into  Eq.  (14.1)  K  becomes  the 
plane-stress  fracture  toughness,  Kc,  ' 

Mira)*  /.\ 

Kc  -  -£-j -  Cl»k)f(|J  •  (14.2) 

Since  the  a/w  ratio  was  identical  for  all  specimens,  the  value 
of  the  finite-width  correction  factor  was  the  same  for  all  tests 
Consequently  the  finite  width  correction  factor  had  no  influence 
on  the  relative  magnitudes  and  trend  of  the  biaxial  load  effect. 

Equation  (14.2)  shows  the  dependence  of  Kr  on  k  as  being 
linear  with  the  factor  (1+k)  increasing  from  one  to  two  as 
k  is  increased  from  zero  to  one.  Therefore  the  value 
would  be  expected  to  approximately  double  over  the  range  of 
k  values  employed  in  these  tests.  The  test  results  shown  in 
Fig.  45  are  in  very  good  agreement  with  this  predicted  trend, 
and  offer  a  further  confirmation  of  the  analytical  portion 
of  this  work. 
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SECTION  IS 


EFFECT  OF  BIAXIAL  LOADS  ON  THE  ANGLE 
OF  INITIAL  CRACK  EXTENSION 

Another  effect  of  load  biaxiality  on  fracture  behavior, 
one  which  is  directly  measurable  and  can  be  compared  with 
the  analytical  predictions,  is  the  influence  of  k  on  the  angle 
of  initial  crack  extension.  The  calculations  of  this  angle, 
based  on  the  maximum  tensile  stress  criterion,  are  displayed 
graphically  in  Fig,  19  for  different  orientations  of  the 
crack  (or  starter  slot).  The  fractured  specimesn  of  all  of 
the  test  series  were  examined  for  comparison  with  the  analyti¬ 
cal  predictions,  and  the  results  are  presented  in  this  section. 

15.1  ANGLE  OF  INITIAL  CRACK  EXTENSION  FOR  THE  ANGLE-CRACKED 

SPECIMENS 

The  effect  of  applied  load  biaxiality  on  Kc  for  the  two 
series  of  tests  on  7075-T6  aluminum  in  which  the  starter  slot 
was  oriented  at  45°  was  discussed  in  the  previous  section. 

The  broken  specimens  from  these  tests  were  examined  to 
determine  the  angle  of  initial  crack  propagation  for  compari¬ 
son  with  the  analytical  predictions  contained  in  Fig.  19. 

As  noted  in  Section  13,  the  material  from  which  these  speci¬ 
mens  were  made  exhibited  a  measurable  amount  of  anisotropy. 

It  was  anticipated  therefore  that  this  could  have  some  in¬ 
fluence  on  the  direction  of  initial  crack  extension.  In 
much  the  same  manner  that  the  average  of  the  two  test  results 
at  each  biaxiality  ratio  minimized  the  effect  of  material 
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anisotropy  on  the  Kc  data  shown  in  Fig.  45,  the  average  of 
the  two  angles  of  initial  crack  extension  for  any  k  value 
should  provide  a  result  relatively  free  of  the  material  ani¬ 
sotropy  effects. 

A  good  appreciation  of  the  influence  of  material  ani- 
sotorpy  on  the  angle  of  initial  crack  extension  can  be  obtain¬ 
ed  from  Fig.  44.  The  analytical  predications  indicate  that, 
for  a«45°  and  k»1.0,  the  crack  angle  should  not  deviate  from 
the  plane  of  the  starter  slot  (0q«Q°).  Figure  44  shows,  how¬ 
ever,  that  the  cracks  propagated  at  an  angle  of  approximately 
±25°  from  the  starter  slot  toward  the  rolling  direction.  If 
the  angles  measured  from  the  two  tests  are  averaged,  they 
approximately  cancel  each  other  leaving  a  net  measured  value 
of  approximately  9o*3°.  This  procedure  was  followed  for  the 
other  three  pairs  of  tests  and  the  results  of  all  eight 
tests  are  summarized  in  TaMe  I.  It  is  seen  from  this  table 
that  the  agreement  between  the  predicted  and  measured  values 
of  6q  is  quite  good,  with  the  only  significant  deviation 
being  for  k»0.5, 

15.2  ANGLE  OF  INITIAL  CRACK  EXTENSION  FOR  HORIZONTALLY 

SLOTTED  707S-T6  SPECIMENS 

Additional  comparison  between  the  predicted  and  measured 
values  of  initial  crack  propagation  can  be  obtained  from  the 
fracture  test  series  discussed  in  Section  13.  For  these 
tests  the  biaxiality  ratio  was  varied  from  zero  to  1.8,  with 
a«90°.  Examination  of  the  predicted  values  of  0Q  (for  a«90°) 


given  in  Fig.  19  shjws  that,  at  the  higher  positive  value? 
o£  k  the  predicted  values  of  9Q  begin  to  increase  rapidly  as 
k  exceeds  values  of  three,  and  approach  values  of  8f-°  to  90° 
as  k  assumes  values  between  five  and  nine.  A  greater  accuracy 
of  prediction  is  difficult  because  of  the  very  high  slopes 
(rates  of  change)  of  the  curves  at  ot-90°.  About  all  that 
can  be  certain  is  that  the  expected  angle  of  initial  crack 
extension  increases  with  increasing  k. 

A  qualitative  confirmation  of  the  above  statement  can 
be  seen  in  Fig.  45,  which  shows  the  effect  of  increasing  k  on 
the  angle  of  initial  crack  extension  for  the  7075-T6  specimens. 
For  these  tests  the  starter  slot  was  perpendicular  to  the 
rolling  direction.  It  is  seen  in  Fig.  45a  (k=0)  that  the 
crack  extended  essentially  parallel  to  the  starter  slot, 
while  in  Fig.  44b  Ck* 1.5)  the  angle  of  initial  crack  extension 
had  increased  significantly.  When  k  was  further  increased 
to  1.8,  it  is  seen  in  Fig.  44c  that  the  crack  very  quickly 
changed  orientation  to  become  more  normal  to  the  direction 
of  maximum  tensile  stress. 

15.3  ANGLE  OF  INITIAL  CRACK  EXTENSION  FOR  PLEXIGLASS  SPECIMENS 

Since  the  material  anisotropy  associated  with  the  7075- 
T6  tests  has  some  effect  on  the  angle  of  initial  crack  ex¬ 
tension  (especially  for  am90°) ,  the  plexiglass  specimens 
were  also  examined  because  they  do  not  exhibit  significant 
anisotropy,  and  because  higher  k  values  could  be  employed  in 
the  tests.  A  large  number  of  tests  were  performed  so  that 
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many  specimens  were  available  for  the  purpose  of  examining  the 
angle  of  initial  crack  extension.  The  results,  however,  for 
any  given  k  value  were  quite  consistent,  consequently,  only 
three  of  the  typical  results  are  presented  in  Fig,  46.  The 
broken  specimen  from  test  B20  (k-0)  is  shown  in  Fig.  46a, 
where  it  is  seen  that  the  crack  extended  with  only  slight 
deviation  from  the  plane  of  the  starter  slot.  Fig, 46b  shews 
the  broken  test  from  test  B17  (k*1.0),  in  which  significant 
deviation  from  the  starter  slot  is  evident.  For  test  B19 

s. 

(k*2.0),  Fig.  46c  shows  that  the  deviation  from  the  plane  of 
the  starter  slot  is  quite  strong,  as  the  crack  attempts  to 
\  align  itself  normal  to  the  direction  of  maximum  tensile  stress. 

All  of  the  plexiglass  specimens  exhibited  similar  tendencies 
|  and  confirm  the  qualitative  agreement  with  the  predicted 

trend  that  was  observed  for  the  7075-T6  specimens  with 
initially  horizontal  crack  orientation.  Additional  observa¬ 
tions  of  the  angle  of  initial  crack  extension  have  been  made 
*  for  PMMA  by  Leevers,  Radon  and  Culver  [42],  who  obtained 

results  in  generally  good  agreement  with  those  presented 
here. 


SECTION  16 


BIAXIAL  FATIGUE  CRACK  GROWTH  RATE  TESTS 


Because  of  the  importance  of  fatigue  crack  growth  laws 
in  determining  the  residual  life  of  structures  and  the  re¬ 
lative  merits  of  many  structural  materials  subjected  to  cyclic 
loading,  one  of  the  objectives  of  this  research  program  was 
to  examine  the  influence  of  biaxial  loads  on  the  fatigue 
crack  growth  rates.  Most  fundamental  studies  of  fatigue 
crack  growth  under  cyclic  loading  conditions  1 4 3 ]  conclude 
that  the  size  of  the  crack-tip  plastic  zone  is  a  critical 
factor  influencing  the  fatigue  crack  growth  rates.  Since 
the  present  study  has  shown  that  increasing  k  values  lead 
to  reduced  crack- tip  plastic  zone  sizes,  it  follows  that 
higher  k  values  should  cause  reduced  fatigue  crack  growth 
rates . 

Experimental  data  about  the  effect  of  load  biaxiality 

on  fatigue  crack  growth  rates  are  currently  quite  limited 

and  somewhat  inconsistent  in  their  assessment  of  biaxial 

effects.  Most  test  results  reported  in  the  literature  were 

performed  at  zero  or  positive  E  values  (IT  •  F.  /F  ) 

min  max 

although  some  exceptions  have  been  reported.  In  tests  con¬ 
ducted  on  a  high  strength  aluminum  alloy,  RR58,  Hopper  and 
Miller  [4*3  reported  a  decrease  in  fatigue  crack  growth  rates 
as  k  increased  from  -1.0  to  1.0.  Kibler  and  Roberts  [35] 
observed  the  same  qualitative  trend  for  6061-T4  and  T6 
aluminum  alloys  as  k  was  increased  from  0  to  0.33.  Leevers, 
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Radon  and  Culver  [4$ ]  reported  fatigue  crack  growth  rates 
in  PMMA  to  be  reduced  by  a  factor  of  two  to  three  as  k  in¬ 
creased  from  0  to  2.0.  For  all  of  the  test  results  just 
ci-ed  the  biaxial  loads  were  in-phase  cyclic  loads. 

On  the  other  hand,  a  variation  of  the  static  lateral 
load  from  k-0  to  k«2.0  appeared  to  have  little  effect  on  the 
fatigue  crack  growth  rates  of  the  more  ductile  PVC  (poly¬ 
vinylchloride)  [45 ) .  Similarly,  liu,  Allison,  Dittmer  and 
Yamane,  observed  negligible  effects  on  the  crack  growth  rates 
of  707S-T73S1  and  2024-T3S1  aluminum  alloys  as  k  was  varied 
from  -1.5  to  1.75  [46],  while  Pook  and  Holmes  found  little 
lateral  load  effect  on  the  growth  rates  of  nickel  alloy 
plates  as  k  was  varied  from  0  to  2.0  [47].  Because  of  the 
lack  of  agreement  between  the  results  of  those  studies  it 
was  considered  desirable  co  perform  some  experiments  to  pro- 
ide  additional  data  about  the  nature  of  biaxial  effects. 

16.1  BIAXIAL  FATIGUE  TESTING  PROGRAM 

In  order  to  gain  initial  insight  into  the  effect  of 
biaxial  cyclic  loads  on  the  fatigue  crack  growth  rates,  a 
preliminary  testing  program  involving  four  specimens  made 
from  thin  sheets  of  2024-T3  aluminum  was  initiated.  For 
these  tests  two  specimens  were  tested  with  the  starter  slot 
parallel  to  the  rolling  direction  (ka0  and  1.0)  and  two  were 
tested  with  the  starter  slot  perpendicular  to  the  rolling 
direction  (k*0  and  1.0).  The  specimens  were  made  identical 
to  those  for  the  static  tests  and  the  loads  were  varied  in 
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a  sinusoidal  manner,  with  both  axis  loaded  in  phase.  The 
load  ratio,  IT,  was  approximately  zero  for  all  four  tests. 

The  test  frequency  was  12  H2  and,  with  the  exception  of  the 
first  test,  the  crack  growth  measurements  were  made  while 
the  test  was  running.  An  arbitrary  initial  cyclic  load 
level  was  imposed  until  the  fatigue  crack  initiated  from  the 
starter  slot  and  the  load  level  was  then  reduced  until  a 
reasonable  crack  growth  rate  was  obtained.  This  load  level 
was  then  maintained  until  the  specimen  fractured.  The  crack 
size  was  recorded  as  a  function  of  the  number  of  cycles  and 
calculations  were  then  made  for  the  crack  growth  rate, 
as  a  function  of  the  applied  stress  intensity  range,  AK. 

The  experimental  results  of  fatigue  crack  growth  rate, 
plotted  on  logarithmic  scale,  versus  AK  is  given  in  Figs. 

48-51,  in  which  each  figure  represents  the  test  results  for 
one  specimen.  Thus  Figs.  48  and  49  represent  the  data  for 
the  two  tests  on  2024-T3  with  the  starter  slot  perpendicular 
to  the  rolling  direction  and  k  equal  to  0  and  1.0  respectively. 
It  can  be  seen  that  there  is  a  significant  biaxial  effect  on 
the  fatigue  crack  growth  rates,  but  determination  of  the 
magnitude  of  the  effect  is  difficult  because  of  the  scatter 
in  the  data.  The  scatter  was  apparently  due  primarily  to  the 
lack  of  experience  in  collecting  fatigue  crack  growth  data 
on  these  biaxial  specimens,  since  it  was  considerably  re¬ 
duced  in  subsequent  tests. 

In  order  to  obtain  a  better  appreciation  for  the  extent 
of  the  biaxial  effect,  a  few  specific  comparisons  will  also 
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be  presented.  In  both  pairs  of  tests,  increasing  the  k  value 
from  zero  to  1.0  caused  a  significant  reduction  in  the 
fatigue  crack  growth  rates.  For  example,  for  the  tests  with 
the  starter  slot  parallel  to  the  rolling  direction,  the  crack 
growth  rates  corresponding  to  the  same  value  of  AK  were 
approximately  1.8x10"^  in/cycle  for  k»0  and  1.4xl0'5  in/cycle 
for  k»1.0.  From  these  results  it  is  seen  that  the  decrease 
in  fatigue  crack  growth  rates  was  approximately  25  percent 
for  an  increase  in  k  from  zero  to  1.0.  The  other  pair  of 
tests  with  the  starter  slot  perpendicular  to  the  rolling 
direction  exhibited  a  greater  dependence  of  the  fatigue 
crack  growth  rates  on  k.  For  example  the  fatigue  crack  growth 
rates  corresponding  to  a  fixed  AK  were  approximately  4.2xl0“5 
in/cycle  ror  k-0  and  2.0x10"^  in/cycle  for  k*1.0,  which  re¬ 
presents  a  change  of  100  percent  in  the  fatigue  crack  growth 
rates.  However,  the  procedure  for  taking  crack  growth  measure¬ 
ments  for  the  k*0  test  was  different  from  the  k«1.0  test  in 
that  the  test  system  was  stopped  while  the  measurements  were 
made.  For  all  subsequent  tests  the  crack  growth  measurements 
were  made  while  the  test  was  running.  A  number  of  subsequent 
tests  were  then  performed  for  the  purpose  of  obtaining  addi¬ 
tional  data  about  the  biaxial  effects. 

These  data  are  presented  in  a  different  way  in  Figs.  52 
d  A 

and  S3,  in  which  ^  is  still  represented  on  the  ordinate 
(vertical  axis)  but  the  half-crack  length,  a,  is  plotted  on 
the  abscissa.  The  data  shown  in  this  manner  for  k«0  and 

1.0  illustrate  quite  clearly  the  biaxial  effect.  The  number 


I 
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of  data  points  have  been  reduced  in  these  figures  by  averaging 
several  adjacent  values,  and  it  is  seen  that  the  scatter 
is  reduced  markedly. 

The  results  of  these  four  tests  are  also  shown  with 
plotted  linearly  on  the  ordinate  and  the  half-crack  length 
plotted  on  the  abscissa.  The  results  are  given  in  Figs.  54 
and  55,  and  again  exhibit  the  significant  biaxial  influence 
already  indicated.  These  results  also  show  very  effectively 
the  acceleration  in  fatigue  crack  growth  rates  that  is 
characteristic  of  these  tests.  Since  these  results  could 
be  easily  misinterpreted  if  different  initial  crack  sizes 
or  different  loads  were  employed  in  the  testing  program, 
results  similar  to  Figs,  54  and  55  are  not  presented  for  the 
subsequent  tests. 

16.2  BIAXIAL  FATIGUE  TESTS  ON  7G75-T6  ALUMINUM 

The  next  test  series  that  was  performed  was  planned  to 
provide  fatigue  crack  growth  rate  clata  over  a  larger  range 
of  k  values.  Thus,  four  additional  tests  were  performed  on 
7075-T6  alumunim  sheets,  0.063  in.  thick,  in  which  the  starter 
slots  were  parallel  to  the  rolling  direction.  The  starter 
slots  were  also  considerably  shorter  than  for  the  prior 
tests  so  that  more  data  could  be  collected.  Also  for  these 
tests,  after  an  initial  set  of  data  were  taken  at  low  crack 
growth  rates,  the  load  was  increased  in  several  steps  so 
that  a  larger  range  of  AK  values  could  be  examined.  The 
cyciic  frequency  for  these  tests  was  decreased  from  15  to 


10  Hz  as  the  test  progressed  and  the  crack  growth  rates  in 
creased.  The  K  value  was  0.1  for  all  four  tests. 


The  effects  of  load  biaxiality  on  the  fatigue  crack 

growth  rates,  as  a  function  of  AK  for  the  four  biaxial!* 

ties,  k«0,  0.5,  1.0  and  1.5  are  shown  in  Figs.  56-S9.  It  is 

seen  from  these  figures  that  the  data  points  exhibited  re* 

asonable  scatter  and  that  a  biaxial  effect  is  clearly  evident. 

As  a  measure  of  the  effect  of  load  biaxiality,  the  fatigue  crack 

growth  ratc-s  for  all  four  tests  at  two  selected  values  of  AK  are  given 

in  Table  II.  The  results  at  both  values  of  AK  exhibited  a  definite 

biaxial  effect,  although  the  primary  effect  seems  to  have 

occurred  between  k«0  and  0.5.  Because  of  the  scatter  that  is 

usually  associated  with  fatigue  crack  growth  rate  testing, 

the  primary  conclusion  from  these  data  is  that,  although 

there  is  a  clear  biaxial  effect,  the  trend  in  the  variation 
da 

of  with  k  will  need  additional  testing  before  it  is  clearly 
understood.  The  results  for  the  first  three  tests,  k*0,  0,5 
and  1,0  have  been  combined  into  one  plot  as  shown  in  Fig.  00, 
where  AK  has  also  been  plotted  on  a  logarithmic  scale.  The 
test  for  k«l.S  was  not  included  because  of  the  very  limited 
amount  of  data  that  had  been  collected  before  the  specimen 
failed  by  breaking  one  of  the  tangs  parallel  to  the  starter 
slot.  The  very  strong  biaxial  effect  is  more  evident  in  this 
figure,  although  the  close  agreement  between  the  results  for 
k*0.5  and  1.0  is  difficult  to  interpret. 

A  different  presentation  of  these  results  is  also  given 
in  Fig.  61.  In  this  figure,  the  half*crack  length,  a,  is 
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plotted  on  the  ordinate  using  a  linear  scale  and  AK  is 
plotted  linearly  on  the  abscissa.  This  presentation  also 
illustrates  clearly  the  magnitude  of  the  biaxial  effect  on 
the  fatigue  crack  growth  rates. 

16.3  BIAXIAL  FATIGUE  TESTS  ON  2024-T3  ALUMINUM 

Since  it  was  considered  likely  that  accurate  fatigue 
crack  growth  rates  are  more  difficult  to  obtain  in  more 
brittle  materials,  three  additional  tests  were  preformed  on 
2024T-3  aluminum  with  biaxiality  ratios  of  0,  0.5  and  1.0. 

The  starter  slots  were  parallel  to  the  rolling  direction  and 
R*0.1  for  all  three  tests.  The  same  general  procedures  estab¬ 
lished  in  the  previous  tests  were  followed  for  these  tests 
also.  The  results  of  the  crack  growth  rates  as  a  function 
of  AK  for  the  three  k  values  are  give  \  in  Figs.  62-64,  where 
these  data  have  been  plotted  on  semi-log  paper.  The  same 
data  were  also  plotted  in  Figs.  65-67  on  a  log-log  plot 
which,  as  was  seen  in  the  7075-T6  results,  tends  to  make  the 
biaxial  effect  less  apparent.  Again  for  direct  comparison 

d  3. 

purposes,  the  values  of  ^  as  a  function  of  k  are  given  in 
the  Table  III.  The  AK  values  of  10  and  40  ksivTn  provided  a 
large  range  of  values  for  additional  comparison  purposes. 

It  is  seen  in  these  results  that  a  significant  biaxial  effect 
exists,  although  the  magnitude  of  the  effect  is  approximately 
30  percent,  rather  than  a  factor  of  10  as  exhibited  for  the 
7075-T6  data.  Also  the  limited  differences  between  the  values 
for  k*0.5  and  1.0  render  interpretation  difficult.  The 
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closeness  of  the  k»0.5  and  1.0  values  for  both  alloys  is  a 
matter  that  can  only  be  a  subject  for  speculation  with  the 
limited  data  that  are  available. 
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PART  III.  CONCLUSIONS 
SECTION  17 

SUMMARY  AND  CONCLUSIONS 


The  plane  problem  of  linear  elasto-statics  is  a  boundary 
value  problem  for  the  solution  of  the  biharmonic  partial 
differential  equation  that  must,  in  addition,  satisfy  given 
conditions  of  stress  or  displacement  that  are  prescribed  on 
the  boundaries  of  the  body.  This  problem  can  be  reformulated 
into  one  of  finding  a  pair  of  analytic  functions  $(z)  and  ft(z) 
of  the  complex  variable,  z,  that  must  satisfy  the  same  pre¬ 
scribed  conditions  of  boundary  stress  or  displacement.  Once 
these  functions  have  been  determined,  the  problem  may  be 
considered  to  be  solved  because  all  of  the  quantities  that 
define  the  mechanical  state  of  the  elastic  body  are  derivable 
from  <t>  and  ft. 

For  the  boundary  value  problem  represented  by  the  biaxially 
loaded  infinite  sheet  with  an  inclined  centrally  located  crack, 
[cf.  Fig.  1],  the  solution  functions  $  and  ft  are  given  by  Eqs. 
(2.8)  and  (2.9).  For  the  horizontally  oriented  crack  they 
are  given  by  Eq.  (2.11),  while  for  the  pair  of  collinear 
cracks,  by  Eqs.  (2.35)  and  (2.37).  In  all  of  these  expressions 
for  $  and  ft,  both  of  the  boundary  load  parameters,  a  and  k, 
appear,  as  they  must,  if  $  and  ft  are  to  represent  solutions  of 
a  boundary  value  problem.  Consequently  it  follows  directly 
from  Eqs.  (2.1)  and  (2.2)  that  the  stress,  displacement,  strain, 
maximum  shear  stress  and  the  strain  energy  density  at  every 
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point  of  the  body  must,  likewise,  show  dependency  on  the 
biaxial  load  parameter,  k.  The  same  must  also  be  true  of 
integrals  of  these  quantities,  yielding  the  elastic  strain 
energy  and  the  elastic  potential  energy  of  the  entire  body. 
When  the  elastic  strain  energy  of  the  body  is  differentiated 
with  respect  to  the  crack  size,  the  resulting  derivative 
representing  the  change  of  the  elastic  strain  energy  with  in¬ 
creasing  crack  size,  must  also  exhibit  a  dependence  on  the 
horizontal  load  parameter,  k,  when  the  body  is  biaxially 
loaded.  Since  every  quantity  defining  the  mechanical  state 
of  the  cracked  body  is  thus  shown  to  be  influenced  by  the 
presence  of  the  horizontal  load,  the  conclusion  that  the 
fracture  behavior  of  the  body  must  also  be  so  influenced  is 
inescapable.  These  conclusions,  which  are  here  drawn  from 
general  considerations  only,  are  clearly  illustrated  by  the 
extensive  analytical  results  appearing  in  the  body  of  the 
report,  and  by  the  experimental  data  that  i,s  presented, 
which  quantitatively  and  qualitatively  confirm  the  findings 
of  the  analysis. 

Stress  and  Displacement  Near  the  Crack-Tip 

The  expressions  for  the  elastic  stress  and  displacement 
in  the  immediate  vicinity  of  the  crack  tip  when  properly 
derived,  will  show  a  dependence  on  the  horizontal  load  para¬ 
meter,  k.  Examination  of  the  series  expansions  of  the 
solution  functions  <t>,  (and  also  for  <t>  ,w)  reveal  that  a 
two-term  approximation  of  the  infinite  series  representations 
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is  necessary  if  proper  account  of  the  presence  of  the  horizon¬ 
tal  load  is  to  be  reflected  in  the  crack-tip  stress  and  dis¬ 
placement  approximations,  [e.g.,  Eqs.  (2.17),  (2.19),  (2.21)- 
(2.23)  for  the  inclined  crack,  Eqs.  (2 . 24) - (2 . 26)  for  the 
horizontal  crack,  and  Eqs.  (2 . 41) - (2 . 47)  for  the  pair  of 
collinear  cracks]. 

When  the  crack  is  horizontal,  the  sec  id  term  of  the 
series  expansion  of  4  is  the  only  term  of  the  expansion  that 
includes  the  parameter  k.  Consequently,  the  usual  one-term 
(singular)  approximation  for  the  crack-tip  stress,  and  also 
for  the  displacement,  is  incapable  of  incorporating  the 
presence  of  the  horizontal  load,  ka,  in  a  biaxial  load  situa¬ 
tion  for  this  crack  geometry. 

Maximum  Shear  Stress  Near  The  Crack  Tip 

The  maximum  shear  stress  about  the  crack  tip  is  shown 
to  be  strongly  influenced  b,  the  presence  of  load  applied 
parallel  to  the  crack  [c.f.  Eqs.  (3. 1) ,  (3. 8) ] .  The  predicted 
pattern  of  variation  of  the  maximum  shear  stress  contours 
with  load  biaxiality,  Fig.  11,  is  confirmed  by  the  photo¬ 
elastic  maximum  shear  stress  isochromatics  that  were  obtained 
for  the  biaxial  load  range  Osksl.O,  shown  by  Fig.  37. 

Because  of  the  relation  between  maximum  shear  stress  and 
plastic  yield,  the  plastic  yield  region  about  the  crack  tip 
should  show  a  qualitatively  similar  pattern  of  variation  with 
load  biaxiality  to  that  shown  by  Fig.  11,  that  is,  as  the 
horizontal  load  varies  from  tension  to  compression  the  extent 


of  the  plastic  yield  region  about  the  crack  tip  should  in¬ 
crease.  (This  predicted  effect  has  been  observed  from  the 
finite-element  (N AST RAN)  elastic-plastic  calculations  con¬ 
ducted  by  Liu  and  Dittmer  [52]  for  a  biaxially  loaded 
cruciform  shaped  specimen.) 

Angle  of  Initial  Crack  Extension 

Utilizing  the  maximum  tensile  stress  criterion,  the 
influence  of  load  biaxiality  on  the  angle  of  initial  crack 
extension  was  calculated  for  different  orientations  of  the 
crack,  with  the  results  summarized  in  Fig.  19.  Test  results 
for  cracks  initially  horizontal  are  shown  in  Fig.  46.  These 
photographs  show  clearly  the  pattern  of  turning  of  the  angle 
of  initial  crack  extension  with  increasing  horizontal  tensile 
load,  in  qualitative  agreement  with  the  predicted  trend  shown 
in  Fig,  19  at  a-Tr/2.  For  the  crack  inclined  at  a*45°,  c.f. 

Fig.  44,  the  measured  angles  of  initial  crack  extension  are  in 
generally  good  agreement  with  the  predicted  values,  rummarized 
in  Table  I. 

Fracture  Load  and  Fracture  Toughness 

A  calculation  of  the  theoretical  fracture  stress,  oc> 
based  upon  Griffith's  crack  instability  hypothesis,  shows 
the  critical  stress  to  be  dependent  on  both  the  load  biaxiality, 
and  the  Poisson  ratio  of  the  material,  c.f.  Eqs.  (9.2)  and 
(9.3),  and  Figures  27,  28  and  29.  In  the  case  of  plane  stress, 
for  materials  with  Poisson  ratio  less  than  one-third,  the 
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breaking  stress  increases  with  increasing  horizontal  tensile 
load.  However  for  values  of  Poisson  ratio  greater  than  one- 
third,  this  trend  is  reversed.  A  similar  patter  of  variation 
holds  for  plane  strain. 

The  results  of  fracture  tests  performed  on  7075-T6  alum¬ 
inum  with  Poisson  ratio  about  0.30,  and  plexiglass  with  Poisson 
ratio  about  0.40-0.45,  are  shown  in  Figures  38  thru  40.  Add¬ 
itional  fracture  test  data  on  PMMA  specimens  (similar  to 
plexiglass)  that  were  obtained  elsewhere  are  shown  in  Fig.  41. 
These  test  results  confirm,  qualitatively,  the  predicted 
influence  of  the  Poisson  ratio  on  the  pattern  of  variation  of 
the  critical  stress  with  variation  of  the  horizontal  load. 

Plane  stress  fracture  toughness  tests  for  determination 

of  K  values,  on  7075-T6  aluminum  sheets  with  the  crack  in- 
c 

dined  at  45°  and  the  horizontal  load  parameter  varying 
between  Osksl.O,  show  a  pronounced  dependence  on  the  load 
biaxialitv.  The  predicted  variation  of  Kc  with  k,  given  by 
Eq .  (14.2),  which  shows  Kc  as  approximately  doubling  in  value 
in  a  linear  manner  as  k  varies  from  zero  to  one,  is  in  ex¬ 
cellent  agreement  with  the  experimental  data  given  in  Fig.  45. 

Fatigue  Crack  Growth  Rates 

The  reduction  of  the  rate  of  fatigue  crack  growth  with 
increase  of  the  horizontal  load,  which  is  implied  by  the  de¬ 
crease  in  size  of  the  crack-tip  plastic  zone  for  such  biaxial 
load  variation,  is  demonstrated  by  the  experimental  data  given 
in  Figures  48  through  67,  and  in  Tables  II  and  III.  These 
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results  are  consistent  with  the  generally  accepted  view  that 
the  size  or  extent  of  the  plastic  yield  region  at  the  crack 
tip  is  a  major  factor  in  influencing  fatigue  crack  growth 
rates.  The  reduction  of  the  crack  growth  rate  as  k  was  varied 
from  0.0  to  0.5  was  observed  to  be  much  greater  than  the  re¬ 
duction  in  rate  that  accompanied  the  change  of  k  from  0.5  to 
1.0.  At  this  time  no  explanation  for  this  surprising  differ¬ 
ence  can  be  offered. 
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SECTION  18 


RECOMMENDATIONS  FOR  FURTHER  RESEARCH 

Although  a  considerable  amount  of  research  effort  has  by 
now  been  directed  toward  examining  the  effects  of  biaxial 
applied  loads  on  the  fracture  characteristics  of  materials, 
many  facets  of  this  effect  have  yet  to  be  clarified  and 
fully  understood.  A  number  of  suggestions  for  further  research 
in  this  area  which  appear  to  be  important  are  indicated  below. 

1.  A  finite-element  analysis  that  would  take  this  particu¬ 
lar  specimen  geometry  into  consideration  should  be  per¬ 
formed.  This  analysis  would  serve  to  provide  more 
accurate  correlations  between  the  anlytical  and  experi¬ 
mental  results.  Such  an  analysis  should  also  lead  to 
further  clarification  of  some  of  the  observed  trends  in 
some  of  the  experimental  results,  such  as  the  interactions 
between  Poisson’s  ratio  and  load  biaxiality  on  the 
fracture  strength. 

2.  The  finite-element  analysis  should  be  extended  to  include 
elastic-plastic  material  behavior  and  subcriHcal  crack 
growth.  An  analysis  of  this  kind  should  provide  a 
greater  appreciation  of  biaxial  load  effects  on  crack- 
tip  plasticity,  and  its  consequent  effect  on  the  fracture 
characteristics.  The  subcritical  crack  growth  capability 
of  the  finite  element  analysis  would  also  serve  to  pro¬ 
vide  more  information  for  comparison  with  the  experi¬ 
mentally-observed  angles  of  initial  crack  extension. 
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The  interaction  between  Poisson's  ratio  and  load  biaxility 

I 

and  their  influence  on  the  fracture  strength  should  be 
examined  further.  The  reversal  of  trend  of  the  fracture  load, 
o  ,  versus  k,  observed  for  aluminum  and  plexiglass  needs  to 
be  examined  in  more  detail  under  both  plane  stress  and 
plane  strain  conditions. 

Fracture  tests  should  be  performed  on  specimens  contain* 
ing  two  collinear  cracks  since  no  tests  of  this  type  have 
been  performed.  It  would  be  desirable  to  perform  tests 
with  horizontal  starter  slots  as  well  as  with  the  starter 
slots  at  an  angle. 

More  fatigue  crack  growth  tests  need  to  be  conducted  to 

del 

examine  in  mo.e  detail  the  biaxial  effects  on  ^  in  the 
range  0<k<1.0,  since  nearly  all  of  the  biaxial  effect 
that  was  observed  in  the  GWU  tests  occurred  as  k  was 
increased  from  zero  to  0.5.  These  test  results  should 
also  be  examined  and  used,  in  conjunction  with  existing 
data,  to  formulate  modifications  to  existing  semi-empirical 
fatigue  crack  growth  rate  laws  that  would  explicitly  in¬ 
clude  k  as  a  factor.  At  present  it  appears  that  the 

d  a 

dependence  on  k  of  jjj  would  be  highly  nonlinear. 

Additional  angle-cracked  tests  should  be  performed  so 
that  crack  orientation  angles  other  than  a»45°  and  90° 
could  be  'tudied.  These  results  would  have  considerable 
practical  significance  since  actual  cracks  in  biaxially- 
loaded  structures  (such  as  pressure  vessels,  etc)  will 
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not  normally  be  aligned  with  the  principal  loading 
directions , 

A  wider  range  of  structural  materials,  including  titanium 
and  steel  alloys,  should  be  tested  in  order  to  determine 
quantitatively  the  relevance  of  the  biaxial  effect  on 
these  structurally  important  materials.  Such  tests 
could  provide  for  greater  variations  in  important 
variables  such  as  crack-tip  plasticity,  subcritical 
crack  growth  and  plane  stress  versus  plane  strain  con¬ 
ditions,  within  the  context  of  a  biaxial  loading  situation. 
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TABLE  I 


Comparison  of  Predicted  and  Measured  Values  for 
the  Angle  of  Initial  Crack  Extension  for  Angle- 
Cracked  Specimens 


BIAXIALITY  RATIO 

0Q (predicted) 1 
(degrees) 

0O (measured) 
(degrees) 

II 

O 

-52 

-48 

k  »  0 . 25 

-46 

-44 

k  *  0.5 

-41 

-27 

k  *  1.0 

-  0 

-  3 

1.  Determined  from  Fig.  19  for  a»45°.  Estimated  accuracy 
±2  degrees  except  for  k=0,2S,  which  was  obtained  by 
interpolation. 

2.  Determined  by  averaging  the  measured  values  of  each 

pair  of  specimens  having  the  same  values  of  k.  .  Estimated 
accuracy  ±3  degrees. 


TABLE  II 


Effect  of  Biaxial  Loads  or.  the  Fatigue  Crack 
Growth  Rates  of  7075-T6  Aluminum 


Biaxiality  Ratio 


§f(in/ cycle) 


AK  *  7ksi/Tn 


AK  *  lOksi  /in 


k  =  0 


3x10 


-5 


9x10 


-5 


k  -  0.5 


3x10 


6.5x10 


k  =  1.0 


3x10 


9x10 


k  =  1.5 


3.5x10 


-6 


3x10 


*5 
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TABLE  III  Effect  of  Biaxial  Loads  on  the  Fatigue  Crack 
Growth  Rates  fo  2024-T3  Aluminum 


1 


Biaxiality  Ratio  ^(in/cycle) 

AK  *  lOksi/in  AK  3  40ksi/in 

k  3  0 

k  3  0.5 

k  3  1.0 


1 . 9x10  ~6  1.3x10 
1 . 5x10  *  ^  1.0x10 
1.4xl0'6  1.1x10 
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4.  Stress  intensity  factor  K7  for  the  inclined 
crack  geometry.  £ 


itr 
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40  SO  60 
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70  80  90 


Fig.  S.  Local  elastic  strain  energy  rate. 


Fig.  11.  Maximum  shear  stress  contours,  horizontal  crack. 


r/a 

FiQ-17,  Error  in  maximum  shaar  stress 
calculation,  shear  panel. 


LOAD  BIAXIALITY  k 

Fig.  29.  Critical  stress  variation  with  load  biaxiality  and  Poisson  ratio 

(plane  stress) 
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Fig.  32.  Illustration  of  fatigue  crack  growth 
rate  versus  AK. 


Figure  33.  Biaxial  specimen  design. 
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Figure  34.  Alternate  biaxial  specimen  design. 
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Fig.  40.  Fracture  toughness  (fracture  load)  for  plastic  sheets 
(George  Washington  University  tests). 


Critical  fracture  load  for  2024-T3  sheets  as  a 
function  of  load  biaxiality.  Crack  perpendicul 
to  rolling  direction. 


(a)  Principal  axis  perpendicular  to  rolling  direction 


(b)  Principal  axis  parallel  to  rolling  direction. 


Fig.  44.  Photographs  of  fractured  angle-crack  specimens 

made  of  707S-T6  aluminum,  0.063"  thick,  k  -  1.0. 
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Fracture  toughness  as  a  function  of  load  biaxiality 
for  angle-cracked  specimens. 
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(b)  k  =  1.5 


(c)  k  =  1.8 


Fracture  paths  of  three  biaxial  test  specimens 
of  7075-T6.  Crack  perpendicular  to  the  roll 
direction. 
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Fig.  49. 
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Fatigue  crack  growth  rate  versus  AK  for  20  24-T3,  k 
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Fig.  60.  Fatigue  crack  growth  rate  versus  AK(logarithmic 
scale)  for  7075-T6.  Crack  parallel  to  rolling 
directici 
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Fig.  64.  Fatigue  crack  growth  rate  versus  AK  for  2024-T3,  k 


CRACK  PARALLEL  TO  ROLLING  DIRECTION 


XjjI 


J _ I _ L  LI  JJLL 


J _ I  I  LUI 1 


10 


100 


AK  (ksi/In) 


Fig.  65.  Fatigue  crack  growth  rate  versus  AK  (logarithmic  scale) 
for  2024-T3,  k  -  0. 
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Fig.  66.  Fatigue  crack  growth  rate  versu: 
for  2024-T3,  k  -  0.5. 
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To  cope  with  the  expanding  technology,  our  society  must 
be  assured  of  a  continuing  supply  of  rigorously  trained 
and  educated  engineers.  The  School  of  Engineering  and 
Applied  Science  is  completely  committed  to  this  ob¬ 
jective. 


